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SUMMARY 


Most  high  temperature  structural  materials  exhibit  a  minimum  in  ductility 
near  one  half  of  their  melting  point.  The  embrittlement  occurs  by  the  * 
nucleation  and  growth  of  cavities  in  the  grain  boundaries.  The  objective  of 
this  investigation  was  to  study  the  phenomenon  of  cavitation  under  conditions  of 
creep-fatigue  interaction  and  under  creep  crack  growth.  The  first  involved 
identifying  and  modeling  the  mechanisms  by  which  cavitation  may  interact  with 
the  classical  stage  I  and  stage  II  crack  growth  in  low  cycle  fatigue.  Creep 
crack  growth  involved  the  growth  of  cavities  in  the  local  stress  and  strain 
field  of  a  crack  tip. 

In  the  creep-fatigue  work,  experiments  were  carried  out  on  model  materials 
to  identify  various  mechanisms  of  failure.  At  least  five  were  found.  The 
dominant  one  depended  on  the  cycle  shape,  the  frequency,  temperature  and  the 
microstructure  of  the  material.  Cavitation  was  found  to  accelerate  the 
initiation  and  propagation  of  fatigue  cracks;  these  were  the  cases  where  creep 
damage  interacted  with  the  cracks  which  are  normally  produced  under  condition  of 
low  cycle  fatigue.  In  other  instances  cyclic  loading  served  to  enhance 
cavitation  damage;  in  this  case  failure  occurred  by  the  gradual  accumulation  of 
cavitation  damage  throughout  the  specimen  rather  than  by  the  initiation  and 
growth  of  a  fatigue  crack.  The  cavitation  mode  of  failure  was  theoretically 
analyzed. 

The  accepted  mechanism  of  creep  crack  growth  is  the  nucleation  and  growth 
of  cavities  in  the  local  stress  field  of  the  crack  tip.  Experiments  were 
carried  out  on  the  nickel  base  superalloy  MERL  76,  to  study  the  influence  of 
loading  rate  on  creep  crack  growth  in  compact  tension  specimens.  We  found  that 
a  slower  loading  rate  decreased  the  fracture  toughness.  The  result  was 
explained  in  terms  of  the  relative  magnitude  of  plastic  flow  and  the  growth  of 
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cavitation  damage  in  front  of  the  crack  tip.  It  was  argued  that  if  the  loading 
conditions  favored  damage  rather  than  plastic  flow,  then  the  toughness  will  be 
lowered.  Local  plastic  strain  and  damage  strain  near  the  crack  tip  were  1 
measured  by  the  distortion  of  a  lattice  of  markers  etched  on  the  specimen 
surface  (using  photo-lithography  techniques).  The  lower  fracture  toughness  was 
found  to  be  associated  with  a  higher  value  of  damage  strain  relative  to  the 
plastic  strain.  A  theoretical  model  was  developed  to  describe  the  results.  The 
model  is  general  and  may  be  used  to  predict  creep  crack  growth  behavior. 

The  grant  provided  whole  or  primary  support  for  one  post-doctoral 
associate,  two  Ph.D.  students  and  one  M.S.  student.  Eight  papers  have  been 
published  in  Metallurgical  Transactions,  Acta  Metallurgica  and  Metal  Science, 
and  one  more  has  been  accepted  for  publication  in  Acta  Metallurgica. 


INTRODUCTION 


Grain  boundaries  have  been  recognized  as  planes  of  weakness  in  structural 
materials  developed  for  high  temperature  applications.  The  embrittlement  occurs 
at  about  the  temperatures  used  in  service.  With  the  development  of  superalloys 
which  have  outstanding  creep  strength,  fracture  rather  than  creep  behavior  has 
become  the  more  important  consideration  in  design. 

The  mechanism  of  failure  at  elevated  temperatures  has  been  recognized  to  be 
the  nucleation  and  growth  of  sub-micron  cavities  in  grain  interfaces.  The 
cavities  prefer  to  form  at  second  phase  particles  such  as  the  carbides  and 
sometimes  the  gamma-prime  precipitates  [1,2].  The  spacing  of  the  cavities, 
therefore,  is  of  the  same  order  as  the  spacing  of  the  precipit?  The  closer 
this  spacing  the  smaller  is  the  ductility  of  the  material.  The  ease  of  cavity 
nucleation  also  plays  an  important  role  in  the  determination  of  time  to  failure. 
It  has  been  hypothesized  that  the  environmental  sensitivity  of  nickel  base 
superalloys  may  be  tied  to  the  nucleation  of  cavities.  One  idea  is  that  oxygen 
diffuses  from  the  surface  along  the  grain  boundaries,  reacts  with  the  carbides 
to  produce  molecules  of  CO  gas  which  promote  the  nucleation  of  cavities  [3], 

The  oxidation  of  the  carbide-matrix  interface  also  weakens  the  boundary  since 
are  generally  less  securely  bonded  to  the  matrix  than  the  carbides  [4], 
e  objective  of  this  investigation  was  to  study  the  mechanisms  of  failure 
under  conditions  of  creep-fatigue  interaction,  and  creep-crack-growth.  The 
approach  was  to  identify  the  mechanisms  and  then  develop  models  for  describing 
failure  in  terms  of  the  loading  variables  and  the  microstructure.  In  the 
following  sections,  the  results  from  each  topic  are  summarized.  Further  details 
of  the  work  are  available  in  publications. 
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In  addition  to  the  work  related  directly  to  creep-fatigue  and  creep-crack 
growth,  additional  papers  on  cavitation  were  also  published  in  the  course  of 
this  investigation,  notably  on  intergranular  fracture  under  aggressive  ' 
environments  and  on  the  description  of  damage  under  multi  axial  loading.  These 
papers  will  not  be  discussed  in  this  report,  but  are  attached  in  the  Appendix. 


CREEP  CRACK  GROWTH 


The  principal  difference  between  stress-rupture,  in  smooth  specimens,  and 
creep-crack-growth,  in  notched  specimens  is  that  in  notched  specimens  the  1 
nucleation  and  growth  of  cavities  is  limited  to  the  crack-tip  region.  In  fact 
the  crack  grows  because  the  cavities  coalesce  ahead  of  the  crack  tip.  A  picture 
of  such  a  mechanism  is  shown  in  Fig.  1.  The  modeling  of  such  a  process  is 
complicated  by  the  following  factors: 

(a)  The  stress  field  near  the  crack-tip  is  non-uniform,  multiaxial,  and 
time  dependent. 

(b)  The  nucleation  and  growth  of  cavities  is  time  dependent  even  under  a 
constant  stress.  Since  the  stress  also  varies  with  time  the 
cavitation  behavior  can  be  analyzed  only  after  making  simplifying 
assumptions. 

(c)  Growth  of  cavities  produces  dilation  in  the  material.  This  can  alter 
the  stress  in  the  crack-tip  region. 

(d)  Since  the  cavities  form  near  the  surface,  environmental  effects  may 
play  an  important  role  in  creep  crack  growth. 

In  the  present  research  program  we  have  completed  three  tasks.  First  a 
simple  model  for  creep  crack  growth  has  been  analyzed.  In  the  model  we  have 
assumed  that  the  cavities  grow  by  diffusion  and  that  the  stress  field  in  front 
of  the  crack  is  determined  by  the  coupling  between  elastic  deformation  and  time 
dependent  cavity  growth.  In  the  second  task  we  have  examined  the  influence  of 
the  rate  of  loading  on  fracture  toughness  in  MERL'76  at  1000K  (1340°F).  We  have 
found  that  fracture  toughness  increases  with  increasing  rate  of  loading.  The 
result  was  explained  in  terms  of  the  relative  magnitude  of  plastic-strain  and 
(cavitation)  damage-strain  near  the  crack-tip.  It  was  argued  that  plastic  flow 
and  cavitation  compete  for  dominance  near  the  crack-tip.  At  slower 


loading-rates,  cavitation  is  dominant  and  the  fracture  toughness  is  low. 

Thirdly,  we  have  developed  a  simple  numerical  model  for  creep  crack  growth  where 
plastic  flow  as  well  as  cavitation  contribute  to  strain  in  the  crack-tip  region. 
The  results  of  the  model  are  in  good  agreement  with  experiments  on  MERL'76.  The 
results  from  each  of  these  tasks  are  now  described. 

Crack  Growth  Model  -  I 

The  growth  of  cavities  in  grain  boundaries  stressed  in  uniform  tension,  has 
been  examined  in  detail  in  the  literature  [5-8].  The  physical  mechanism  of 
cavity  growth  is  illustrated  in  Fig.  2.  Cavities  grow  by  the  transfer  of  atoms 
from  cavity  surfaces  to  the  grain  interface.  As  a  result  the  cavities  grow  in 
size  and  the  matter  "plated"  onto  the  grain  boundary  is  accomodated  by  rigid 
body  displacement  of  the  two  half  crystals.  If  the  cavities  are  assumed  to  have 
an  equilibrium  shape,  then  the  grain  boundary  damage,  Ar,  defined  as  the 
separated  fraction  of  the  boundary  area,  the  normal  displacement  of  the 
crystals,  u,  and  the  average  cavity  spacing,  X,  are  related  to  each  other  by 
geometry.  The  detailed  expressions  are  given  in  Ref.  9.  If  we  normalize  u  with 
respect  to  X,  so  that: 


u 


(1) 


then  the  rate  of  growth  of  cavities  (expressed  in  terms  of  u)  can  be  related  to 
the  (normalized)  applied  stress,  aQ  by: 

du 


—  =  (on—  a*) 
d?  0 


(2) 


In  (2),  d* 


is  the  threshold  stress  for  the  nucleation  of  cavities  and  t  is  a 


normalized  unit  of  time  which  depends  on  temperature  through  the  boundary 
diffusion  coefficient.  The  stress  in  equation  (2)  has  been  normalized  with 
respect  to  the  Young's  Modulus.  Equation  (2)  is  the  simple  rather  than  the 
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exact  version  of  the  growth  equation.  The  latter  is  discussed  in  Ref.  9. 

Recalling  that  the  damage,  is  related  to  u  through  geometry,  the 
fracture  conditon,  defined  as  the  point  where  neighbor  cavities  grow  to  touch 
each  other,  can  be  defined  in  terms  of  u.  It  may  be  shown  that  this  critical 
displacement,  uc  =  0.4. 

The  stress  field  near  a  crack-tip  is  non-uniform  and  the  displacement 
required  to  accomodate  cavity  growth  will  generate  elastic  back  stresses,  which 
will  tend  to  relax  the  stress  field  at  the  crack-tip.  That  is  the  problem  which 
we  have  analyzed  in  detail.  One  set  of  results  are  shown  in  Figure  3.  Growth 
of  cavities,  which  is  characterized  by  the  displacement  u,  relaxes  the  stress 
field.  When  u  =  uc  at  the  crack-tip,  the  crack  advances  by  one  unit  spacing. 

The  calculation  leads  to  the  following  interesting  results: 

(i)  Crack  velocity  reaches  a  steady  state  after  a  readjustment  of  the 
stress  field  near  the  crack-tip. 

(ii)  For  a  crack  moving  in  the  steady  state,  the  stress  field  remains 
constant. 

e\ 

(iii)  The  steady  state  crack-velocity,  a,  is  proportional  to  Kj,  where  Kj  is 
the  applied  stress-intensity  factor. 

(iv)  A  characteristic  damage  zone,  of  size  Z^,  develops  in  front  of  the 
crack-tip.  also  depends  on  Kj. 

(v)  Crack  will  not  propagate  if  Kj<Kt^.  The  threshold  is  related  to  the 
threshold  stress  for  cavity  nucleation. 

The  equations  for  the  above  properties  of  the  model  are  given  below: 


K.h  =  (0.4  AEo*)1/2 


'V  A  /.W  A  V.  AV.  A  A  -• 


•V*  V  V 


1 
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a  =  7xlCr  — - £  when  KT >K. , 

EkT  1?  1  th 


a  =  0 


when  KjSKth 


In  equations  (3)  and  (4),  E  is  the  Young's  modulus  and  Q  is  the  atomic  volume. 

2 

Although  the  equations  agree  quite  with  experiment  in  threshold  behavior  and  Kj 
dependence  of  crack-velocity,  the  measured  crack  velocit  es  [10,11]  appear  to  be 
about  two  to  three  orders  of  magnitude  slower  than  those  calculated.  The 
discrepancy  may  arise  from  (i)  the  assumption  that  the  crack  moves  in  a  flat 
plane  rather  than  in  a  zig-zag  manner  as  would  be  the  case  for  an  intergranular 
path  in  a  polycrystal  and,  (ii)  the  influence  of  plastic  flow,  which  has  been 
disregarded  in  the  model,  on  the  crack-tip  stress  field.  A  comparison  between 
theory  and  experiment  is  given  in  Ref.  9. 

Experiments  to  Measure  Fracture  Toughness  of  MERL'76  at  1000K 


We  have  measured  fracture  toughness  as  the  work  of  fracture  required  to 
initiate  and  propagate  an  intergranular  crack  from  a  notch  in  compact-tension- 
specimens.  The  concept  underlying  the  experiments  is  illustrated  in  Fig.  4(a) 
and  (b).  Deformation  near  the  crack-tip  has  two  components:  plastic  flow,  and 
grain  boundary  cavitation  damage.  Plasticity  concentrates  along  the  shear  lobes 
where  the  shear  "tress  is  maximum.  Cavitation,  on  the  other  hand  would  localize 
immediate  ahead  of  the  crack-tip  since  it  responds  to  the  maximum  tensile 
stress.  Whe.oas  plastic  flow  will  blunt  the  crack-tip,  cavitation  will  move  the 
crack  forward  and  keep  it  sharp.  Our  view  was  that  the  relative  magnitude  of 
plastic  strain  and  damage  strain  would  influence  the  measured  fracture  toughness 
of  the  material.  Their  relative  magnitude  could  be  changed  through  the  loading 
rate.  We  expected  that  the  slower  loading  rates  would  enhance  the  relative 
significance  of  damage  and,  therefore,  lower  the  fracture  toughness. 


The  experimental  results  shown  in  Fig.  4(d)  are  in  agreement  with  the  above 
hypothesis.  In  order  to  confirm  the  hypothesis,  the  damage  strain  and  the 
plastic  strain  were  measured  directly.  The  technique  was  to  etch  markers  on  the 
surface  at  a  spacing  of  6ym.  Plastic  strain  was  defined  as  the  deformation  of 
the  grains,  and  was  measured  as  the  distortion  of  the  etched  grid  on  the  scale 
of  one  grain  size.  Cavitation  strain  was  defined,  and  measured,  as  the  average 
separation  of  the  grain  interfaces  per  unit  distance.  A  picture  of  a  situa  n 

where  the  damage  strain  is  significant  but  the  plastic  strain  is  negligible  - 
shown  in  Fig.  4(c).  A  plot  of  the  ratio  of  the  plastic  strain  to  damage  st  .1 
is  given  in  Fig.  4(d);  note  that  the  increase  in  the  fracture  toughness  is 
directly  correlated  to  this  ratio. 

An  approximate  computer  analysis  was  developed  to  describe  the  results  by  a 
quantitative  model.  The  model  has  been  described  quantitatively  in  Ref. 12.  The 
details  of  the  model  will  not  be  discussed  here.  The  results  of  the  model  were 
in  good  agreement  with  experiment  as  illustrated  in  Fig.  4(d). 


CREEP-FATIGUE  INTERACTION 


Experiments 

Creep-fracture  is  characterized  by  the  nucleation  and  growth  of  cavities  at 
grain  boundaries,  while  low-cycle-fatigue  is  characterized  by  the  Stage  I  crack 
initiation  and  Stage  II  propagation  of  a  transgranular  crack.  Separately  both 
failure  processes  are  well  understood  but  the  interactions  between  them  appear 
to  be  complex.  A  linear  damage  rule  which  assumes  that  the  fatigue  portion 
(defined  by  the  number  of  cycles)  and  creep  portion  (defined  as  the  time  spent 
at  maximum  stress)  add  to  cause  failure  has  proved  to  be  too  simplistic.  The 
data  when  plotted  by  this  rule  show  a  very  large  scatter  suggesting  that  the 
interaction  mechanisms  are  non-linear  and  that  more  than  one  type  of 
interactions  may  be  possible. 

Our  approach  to  the  above  dilemma  was  to  carry  out  systematic  experiments 
to  characterize  the  mechanisms  of  creep-fatigue  interactions.  Aluminum-five 
percent  magnesium  alloy  was  chosen  for  the  study  because  it  has  been  shown  to  be 
very  resistant  to  intergranular  creep-fracture  under  monotonic  loading. 
Experiments  were  done  on  smooth  specimens  with  a  uniform  gage  section  (as 
opposed  to  hour-glass  shaped  specimens)  and  in  ultra-high-vacuum  to  reduce  the 
effects  of  the  environment  to  a  minimum.  The  temperature  and  the  cycle  shape 
were  varied,  but  the  plastic  strain  size  and  the  frequency  were  held  constant. 

It  was  felt  that  frequency  and  temperature  are  complementary  parameters  which 
may  be  normalized  with  respect  to  each  other. 

The  results  of  the  experiments  are  shown  in  Fig.  5b.  The  cycle  shapes 
which  were  used  are  drawn  in  Fig.  5a.  Metallography  and  fractography  revealed 
at  least  five  different  failure  mechanisms  of  creep  fatigue  interaction.  They 
are  designated  by  (a)  thru  (e)  in  Fig.  5b.  Each  one  is  described  below. 

Further  details  and  the  micrographs  are  contained  in  Ref.  13. 


(a)  This  was  the  classical  low-cycle-fatigue  mode  of  failure  which  is  produced 
by  the  transgranular  initiation  of  a  Stage  I  crack  and  propagation  of  a 
Stage  II  transgranular  crack.  The  propagation  is  characterized  by 
striations  on  the  fracture  surface.  The  spacing  of  the  striations  is 
believed  to  be  related  to  the  crack-tip-opening-displacement  [14]. 

(b)  This  failure  mode  occurred  for  all  three  cycle  shapes  at  423K,  and  at  473K 
for  the  equal  ramp  cycle  shape.  It  differed  from  (a)  in  that  crack- 
initiation  changed  transgranular  to  intergranular.  But  propagation 
remained  transgranular  and  was  marked  by  striations. 

(c)  This  failure  mode  dominated  at  higher  temperatures  when  the  cycle  shape 
contained  a  tensile  hold-time.  In  this  case  the  failure  mode  was  totally 
different.  Instead  of  crack  initiation  and  propagation,  intergranular 
creep  damage  was  spread  throughout  the  specimen.  Interrupted  tests  under 
the  same  conditions  revealed  that  the  damage  accumulated  gradually  with 
fatigue  cycles.  The  damage  was  of  1 r'  type  of  cavities  as  shown 
schematically  in  Fig.  6. 

(d)  In  this  case,  as  in  (c),  failure  was  caused  by  cumulative  creep  damage,  the 
difference  being  that  'wedge'  type  of  grain  boundary  cracks  were  formed. 

The  difference  between  'r'  and  'wedge'  cavities  is  illustrated  in  Fig.  6. 
The  cumulative  nature  of  the  damage,  which  was  spread  throughout  the 
specimen,  was  confirmed  through  interrupted  tests. 

(e)  Here  a  macroscopic  instability  and  flow  localization  developed  in  the  gage 
length  of  the  specimen.  It  is  possible  that  the  instability  was  caused  by 
microstructural  degradation  of  magnesium-aluminide  precipitates  at  573K. 

Theory 

The  cumulative  damage  mechanism  of  creep-fatigue  interaction  for  the 

'wedge*  type  of  intergranular  cavities  was  formulated  and  analyzed  in  detail. 


This  is  the  type  (d)  failure  described  in  the  previous  seciton.  The  physical 
basis  of  the  analysis  and  the  cogent  results  are  described  here.  The  details 
are  given  in  Ref.  15.  * 

The  crucial  role  of  grain  boundary  sliding  in  the  growth  of  wedge  cracks 
may  be  appreciated  from  the  schematic  in  Fig.  6(a).  The  basis  of  the  model  then 
is  to  correlate  the  applied  strain  to  the  local  sliding  at  grain  boundaries,  and 
the  sliding  to  the  wedge  crack  damage.  First,  wedge  damage,  Aw,  is  defined  by 
the  equation: 


w 


(5) 


where  w  is  the  length  of  the  wedge  crack  and  L  is  the  grain  size.  The  sliding 
displacement,  v,  at  the  triple  junction  where  the  wedge  crack  initiates  is  also 
normalized  with  respect  to  the  grain  size,  so  that  at: 

v  .  "sliding  (6) 


Then  a  simplifying  assumption  is  made  that  A  and  v  are  linearly  related  through 

w 

a  material  parameter,  a: 


a 


(7) 


The  connection  between  the  applied  strain  eg  and  the  sliding  strain,  v,  is 
complex  but  has  been  extensively  studied  experimentally  and  theoretically. 
Defining  a  quality  6  such  that: 


(8) 


it  may  be  shown  that  B  varies  sigmoidally  with  strain-rate  as  shown  in  Fig.  7. 
At  slow  strain  rates  8  approaches  0.5.  As  the  strain-rate  increases  there  is  a 
transition  until  sliding  becomes  zero  at  very  high  strain-rate.  The  important 


is  a 


.  * 


parameter  is  the  strain-rate  at  which  the  transition  is  felt.  It,  ew  , 


material  parameter  and  has  been  found  to  depend  on  the  microstructure  of  the 
grain  boundary  by  the  following  equation:  ‘ 


*  0.27YH  6Dfa 

ew  =  ,V_  7727 


(9) 


kT  fbp  L 

Here  Y  is  the  yield  stress,  p  is  the  size  of  the  particles  in  the  boundary  and 
fb  is  the  area  fraction  of  the  particles  in  the  boundary.  In  cases  where  these 
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parameters  cannot  be  characterized  we  have  shown  that  ew  may  be  experimentally 


measured  with  the  use  of  internal  friction  measurements  of  grain  boundary 
sliding  rate.  These  details  are  discused  in  Refs.  13  and  17. 

Equations  (7)  and  (8),  when  combined  give  the  result  that: 

dA. 


de. 


=  a 


(10) 


where  B  =  B(ce/ew*). 


For  a  slow-fast  cycle  shape,  the  value  of  8  in  the  tensile  half  of  the 
cycle  will  be  greater  than  in  the  compressive  half  of  the  cycle,  because  the 


tensile  strain-rate  is  slower.  If  the  tension  strain-rate  is  and  the 


compression  strain-rate  is  ec,  then  the  damage  per  cycle  will  be  given  by: 


=  a6('t)  '  ^c* 


(ID 


The  first  term  in  (11)  is  the  positive  damage  caused  by  the  tensile  portion  of 
the  cycle.  The  second  term  is  negative  because  compression  would  tend  to  heal 
the  damage  produced  in  tension.  The  coefficient  y  has  been  asserted  to  account 
for  the  hysteresis  in  the  recovery  of  damage.  The  cycles  to  failure,  Nf=l/A^, 
which  leads  to  the  result: 


Ae_Nx  = 


1 


P  f  a(8(et)  -  3( ee ) 1 


(1?) 


1 


The  result  in  (12)  shows  that  N^.  «  l/Ae^  when  cumulative  creep  damage  is  the 
mechanism  of  failure  (in  low-cycle-fatigue,  Coffin-Manson  equation  predicts  a 
parabolic  power- law-dependence).  Equation  (12)  does  contain  two  adjustable 
parameters  a  and  y  which  can  be  fixed  by  simple  experiments. 

The  application  of  the  above  result  to  our  own  measurements  on  aluminum  and 
to  published  data  for  austentic  stainless  steels  has  been  discjssed  in  Ref.  15. 
The  agreement  is  quite  good.  The  application  of  the  analysis  to  nickel  base 
superalloys  has  been  less  successful  suggesting  that  cumulative  creep  aamage  may 
not  be  the  primary  failure  mode  in  superalloys  at  temperatures  encountered  by 
disks  in  gas-turbine  engines.  We  believe  that  the  intergranular  crack 
initiation  mechanism  (described  by  (b)  in  Fig.  5b)  is  more  relevant  in 
nickel-base  superalloys.  The  sensitivity  of  superalloys  to  oxygen  environments 
could  also  complicate  the  prediction  of  life  by  the  simple  analysis  as  described 
in  this  report. 
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FIGURE  CAPTIONS 


Cavities  forming  ahead  of  a  crack-tip  of  an  advancing  creep  crack  in 
MERL'76  at  1000K. 

The  diffusional  transport  mechanism  of  cavity  growth.  Atoms  are  removed 
from  the  cavity  surface  and  inserted  into  the  grain  boundary  adjacent  to 
the  cavities. 

Numerical  results  from  Ref.  9,  showing  the  relaxation  of  stress  just  ahead 
of  the  crack  tip  with  the  growth  of  cavities.  U  is  a  measure  of  cavity 
growth.  Eventually  the  crack  moves  at  a  steady  state. 

Model,  theory  and  experiment  for  the  fracture  toughness  of  MERL'76  at 
1000K.  Details  given  on  page  8. 

(a)  The  cycle  shapes  employed  in  experiments  to  identify  the  mechanisms  of 
creep-fatigue-interaction,  (b)  The  influence  of  cycle-shape  and 
temperature  on  fatigue  life  and  on  the  failure  mechanism;  five  failure 
mechanisms  were  identified.  Taken  from  Ref.  13. 

A  schematic  of  two  types  of  creep  damage. 

Variation  of  the  ratio,  8,  of  strain  due  to  sliding  to  the  applied  strain 
as  a  function  of  the  applied  strain-rate  eg  [16]. 
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The  diffusional  transport  mechanism  of  cavity  growth.  Atoms  are  removed 
from  the  cavity  surface  and  inserted  into  the  grain  boundary  adjacent  to 
the  cavities. 
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Numerical  results  from  Ref.  9,  showing  the  relaxation  of  stress  just  ahead 
of  the  crack  tip  with  the  growth  of  cavities.  U  is  a  measure  of  cavity 
growth.  Eventually  the  crack  moves  at  a  steady  state. 


MODELING  OF  THE  CREEP  CRACK  GROWTH  PROCESS 


(a)  HIGH  TOUGHNESS 


(b)  LOW  TOUGHNESS 


»  Direction  of 
Crock  Propagation 


c 

E 

w 


■5  Li - 1 _ 1 - 1 - i_l 

S.  0.01  0.1  1  10  100 

Loading  Rate  (MPa  ^m/s) 


(C)  DAMAGE  DOMINATED  DEFORMATION  (d)  INFLUENCE  OF  LOADING  RATE  ON  DAMAGE 


4.  Model,  theory  and  experiment  for  the  fracture  toughness  of  MERL'76  at 
1000K.  Details  given  on  page8. 
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(a)  The  cycle  shapes  employed  in  experiments  to  identify  the  mechanisms  of 
creep-fatigue-interaction,  (b)  The  influence  of  cycle-shape  and 
temperature  on  fatigue  life  and  on  the  failure  mechanism;  five  failure 
mechanisms  were  identified.  Taken  from  Ref.  13. 


Variation  of  the  ratio,  S,  of  strain  due  ,to  sliding  to  the  applied  stra 
as  a  function  of  the  applied  strain  rate  e  [16]. 
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APPENDIX 


The  following  two  papers  discuss  the  environmental  effects  on  cavitation 
and  cavitation  under  multiaxial  loading.  They  are  included  here  because  they 
have  not  been  discussed  in  the  text.  Both  papers  are  particularly  relevant  to 
nickel  base  alloys  because  (a)  they  are  known  to  be  environment  sensitive  and 
(b)  they  exhibit  creep  crack  growth  which  involves  the  growth  of  cavities  in  a 
multiaxial  stress  field  ahead  of  the  crack  tip. 
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Abstract — Intergranular  cavitation  produces  dilation  in  the  specimen.  It  can  also  enhance  the  rate  of 
creep  as  demonstrated  convincingly  by  experiments  of  Needham  and  Greenwood  [I  J  Metals  Sci.  9.  258 
(1975).  Expressions  for  cavitation  enhanced  creep,  and  for  cavitation  induced  dilation  arc  derived.  The 
latter  depends  only  on  geometry,  but  the  first  is  a  consequence  of  the  constrained  cavity  growth 
mechanism.  That  mechanism  involves  redistribution  of  stress  from  regions  which  have  cavities  to  those 
which  do  not;  as  a  result  the  stress  on  the  constraining  regions  is  increased  leading  to  an  increase  in  the 
measured  creep-rate.  A  comparison  of  compressive  vs  tensile  tests  is  suggested  as  a  means  of  dis¬ 
tinguishing  the  enhancement  in  creep-rate. 

Resume — La  cavitation  intergranulaire  produit  des  dilatations  dans  I'echantillon.  Elio  peut  egalement 
augmenter  la  vitesse  de  fluage  comme  l’ont  demontre  dc  maniere  convaincante  les  experiences  de 
Needham  et  Greenwood  [I]  Metals  Sci.  9,  258  (1975).  Nous  avons  obtenu  des  expressions  pour  le  lluage 
induit  par  cavitation  et  pour  la  dilatation  induite  par  cavitation.  Cette  derniere  ne  depend  que  de  la 
geometrie.  mais  la  premiere  est  une  consequence  du  mecanisme  de  croissance  forcec  des  cavites.  C'c 
mecanisme  implique  une  redistribution  de  la  contrainte  dcs  regions  presenlant  des  cavites  vers  cellos  qui 
nen  ont  pas;  il  en  resulte  une  augmentation  de  la  contrainte  dans  les  regions  contraintes  provoquant 
une  augmentation  de  la  vitesse  de  fluage  expcrimentale.  Nous  suggerons  une  comparaison  des  essais  de 
compression  et  de  traction  afin  de  mettre  en  evidence  Taugmenlation  de  la  vitesse  de  linage. 

Zusammenfassung  -Inlergranutare  Hohlraumbildung  fiihrt  zu  eincr  Aufweitung  der  Probe.  Sie  kann 
also  die  Kriechgeschwindigkeil  crhbhen.  wie  die  Experimente  von  Needham  und  Greenwood  [1] 
Metals  Sci.  9,  258  (19751  iiberzeugend  beweisen.  Das  hohlraumbcschleumgte  Kriechen  und  die 
hohlraumbedingte  Aufweitung  werden  mathematisch  beschricben.  Die  Aufweitung  hangt  nur  von  der 
Geometrie  ab.  das  hohlraumbeschieunigte  Kriechen  jedoch  ist  eine  Folge  des  eingcschranktcn  Hohl- 
raumwachstums.  Bei  diesem  Mcchanismus  werden  Spannungen  aus  Bereichcn  mil  Hohlraumen  in  Ber- 
eiche  ohne  Hohlriiume  verlagert;  als  Folge  wird  die  Spannung  auf  die  einschriinkenden  Bcreiche  ver- 
grbBert,  welches  zu  einem  Anstieg  in  der  gemessenen  Kriechgeschwindigkeil  fiihrt.  Fiir  die  Ermittlung 
des  Ansticgs  in  der  Kriechgeschwindigkeil  wird  vorgeschlagen.  Druck-  und  Zugversuchc  miteinander  zu 
vergleichen. 


INTRODUCTION 

Deformation  of  a  cavitating  solid  can  be  described  by 
two  scalar  quantities:  an  effective  strain  ee.  which  is  a 
measure  of  the  change  in  shape  without  a  change  in 
volume,  and  a  cavitation  strain  e0,  which  is  the  frac¬ 
tional  increase  in  volume  due  to  internal  cavities. 
Intuition  suggests  that  in  multiaxial  loading  the  mean 
stress,  or  pressure,  should  influence  ea  but  not  er,  yet 
Needham  and  Greenwood  [I],  and  Dyson  et  ul.  [2] 
have  found  that  the  effective  strain-rate  does  not 
depend  merely  on  the  effective  stress,  as  found  in  plas¬ 
tic  flow  of  metals,  but  also  on  the  state  of  the  applied 
stress.  They  found  that  when  creep  tests  were  done  on 
copper  in  a  variable  hydrostatic  pressure,  then  the 
rate  of  creep  was  influenced  by  the  pressure,  even 
when  the  effective  stress  was  kept  constant.  The  press¬ 
ures  were  too  small  to  influence  the  self  diffusion  coef¬ 
ficient.  Furthermore,  the  change  in  the  creep-rate  with 
pressure  appeared  to  be  unrelated  to  the  extent  of 
cavitation  (as  measured  by  the  change  in  density). 


Dyson  et  ai.  in  their  work  with  stainless  steel,  found  a 
difference  in  the  creep-rate  between  a  tension  test  and 
a  torsion  test,  for  the  same  effective  stress;  but  the 
difference  appeared  only  when  cavitation  was  present. 
In  combination,  the  two  sets  of  results  imply  that 
cavitation  is  a  necessary  condition  for  this  effect, 
although  the  magnitude  of  the  effect  is  not  related  to 
the  extent  of  cavitation.  This  means  the  constitutive 
equation  for  cavitation  enhanced  creep  should  depend 
on  the  stress-state  but  may  not  necessarily  include  a 
state  variable  related  to  cavitation,  such  a.std. 

Dyson  et  at.  proposed  that  the  above  phenomenon 
can  be  explained  in  terms  of  the  constrained  cavity 
growth  mechanism  [3  5],  Their  concept  is  illustrated 
in  Fig.  1.  If  cav  ities  form  preferentially  on  grain  facets 
aligned  nearly  perpendicular  to  the  maximum  positive 
principal  stress  and  if  the  cavities  grow  by  diffusion 
then  grains  such  as  A  and  B  must  move  apart  from 
each  other  to  accommodate  the  volume  of  material 
which  is  displaced  by  cavity  growth.  They  will  be  con¬ 
strained  from  doing  so  if  the  creep-rate  in  the  neigh- 
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Fig.  1.  A  schematic  illustrating  the  constrained  growth  of 
cavities  in  a  biaxially  stressed  specimen.  The  dilation  due 
to  cavity  growth  on  transverse  boundaries  must  be  acco¬ 
modated  by  plastic  flow  in  neighbouring  grain  matrix. 


boring  grains  is  slower  than  the  displacement  rate 
required  to  accommodate  cavity  growth.  The  dis¬ 
placement  rate  then  becomes  limited  by  the  rate  of 
creep  in  adjacent  regions,  such  as  M.  (It  must  be 
remembered  that  there  are  different  regimes  of  behav¬ 
ior  depending  upon  the  stress,  temperature  and 
microstructure.  In  general,  low  stress,  high  tempera¬ 
ture,  small  cavity  spacing,  and  a  long  time-to-failure 
favors  the  constrained  growth  mechanism  [3-5].)  This 
mechanism  gives  a  qualitative  explanation  of  why 
cavitation,  but  not  the  size  of  the  cavities,  may 
influence  the  creep-rate.  In  Fig.  1,  for  example,  the 
load  on  boundary  AB  would  be  transferred  to  region 
M  leading  to  an  increase  in  the  rate  of  creep.  Yet,  the 
magnitude  of  the  effect  will  not  depend  on  cavity  size. 
The  only  condition  for  load  transfer  is  that  if  an  equal 
tensile  stress  were  to  be  applied  to  region  AB  and 
region  M,  isolated  from  each  other,  then  AB  would 
elongate  at  a  much  faster-rate  than  M.  It  has  been 
shown  [6,7]  that  the  displacement  rate  of  AB 
depends  more  strongly  on  cavity  spacing  than  on 
cavity  size. 

When  cavities  form,  then  one  of  the  conditions  for 
pure  plastic  flow,  that  deformation  should  occur  at 
constant  volume,  is  no  longer  satisfied.  In  long  term 
fracture  experiments,  done  under  conditions  which 
approach  the  service  conditions,  the  volumetric  strain, 
in  fact,  can  become  the  dominant  strain  [8, 9],  Thus  it 
becomes  necessary  to  distinguish  between  the  shape- 
change  component  and  the  volumetric  component  of 
the  total  strain,  before  the  usual  multiaxial  stress/ 
strain-rate  equations  for  power-law-creep  [2]  can  be 
applied. 

This  paper  consists  of  four  parts.  First  a  procedure 
for  separating  the  creep  or  the  effective  strain  and  the 
cavitation  strain  is  presented.  Next,  equations  are  de¬ 


veloped  for  describing  the  cavitation-strain  in  terms 
of  the  geometry  and  the  distribution  of  cavities  in  the 
grain  boundaries.  Third,  a  power-law  constitutive 
equation,  in  the  presence  of  growing  cavities,  is  de¬ 
scribed.  Finally,  a  possible  form  for  the  Monkman 
Grant  equation  for  use  under  multiaxial  straining  is 
given.  The  second  and  the  third  items  are  not  new  in 
concept.  Harris  et  al.  [8]  have  also  considered  the 
creep  strain  arising  from  the  growth  of  cavities. 
Although  the  basic  premise  of  their,  and  our 
approach  is  similar,  there  are  important  differences. 
Their  formulation  begins  from  the  equations  by  Hull 
and  Rimmer  [10]  and  Speight  and  Harris  [1 1]  for  dif- 
fusional  cavity  growth.  It  now  appears  that  the  con¬ 
strained  cavity  growth  mechanism  is  more  relevant  to 
long-term  creep  fracture  experiments.  In  this  mechan¬ 
ism,  the  rate  of  cavity  growth  is  related  to  tensile 
strain,  rather  than  to  time  as  dictated  in  the  Hull  and 
Rimmer  approach.  Also,  in  the  present  work  the 
dependence  of  the  rate  growth  of  cavitation  strain  on 
the  following  parameters:  cavity,  spacing,  grain-size 
and  the  di-hedral  angle,  is  explicitly  defined. 

The  third  problem,  which  pertains  to  the  prescrip¬ 
tion  of  a  constitutive  equation  for  cavitation 
enhanced  effective  strain-rate  has  recently  been  con¬ 
sidered  by  Dyson  et  al.  [2],  The  major  difference 
between  their  results  and  the  equations  derived  here 
lies  in  theit  assumption  that  deformation  occurs  at 
constant  volume;  the  use  of  equation  (2)  in  their 
paper  is  contingent  upon  this  assumption.  While  the 
assumption  may  be  reasonably  valid  when  fracture 
occurs  at  large  ductilities,  it  is  misleading  in  long  term 
experiments  when  the  cavitation  strain  can  become  a 
significant  fraction  of  the  total  strain  [9].  The  analysis 
given  here  distinguishes  between  creep  strain,  which  is 
defined  for  deformation  at  constant  volume,  and  the 
cavitation  strain. 

THE  STRAIN-STATE 

Any  strain  state  can  be  described  by  the  three  prin¬ 
cipal  strains.  « e2  and  e3.  and  by  the  orientation  of 
the  axes  along  which  these  strains  are  prescribed  In 
general  then 

’« ,0  o' 

[«]=  0  0  (1, 
_<)  0  fj 

The  strain  due  to  change  in  volume.  and  that  due 
to  change  in  shape  r,„  can  be  separated  as  follows 

f,  =  f,  +<;  I  Jl 

and 

«,  =  ?\  i  M*l  —  « 2  I  ’  *  0;  -  *  il; '  •  I*  ,  <  ,  I  I  1 1 

The  strain  is  related  to  the  change  m  densiiv  h\ 
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3  cos  0  +  cos3  (I). 


(SI 


where  the  approximate  form  holds  for  small  volu¬ 
metric  strain.  Here  A (>,p  is  the  fractional  change  in 
density. 

A  comparison  of  equations  (2)  and  (3)  shows  that  e, 
and  <„  are  independent  of  each  other  and,  therefore, 
measure  different  effects  of  deformation.  Tnis  also 
means  that  if  both  (a  and  «,  are  non-zero  then  at  least 
two  principal  strains  must  be  measured  in  order  to 
distinguish  between  them.  As  an  illustration,  consider 
a  uniaxial  test.  Here  we  must  measure  either  the 
transverse  strain,  e„  or  the  change  in  density,  in  ad¬ 
dition  to  the  axial  strain,  e.,  in  order  to  distinguish 
between  ea  and  €,.  Recognizing  that  two  of  the  princi¬ 
pal  strains  in  a  uniaxial  test  are  equal,  i.e.  the  trans¬ 
verse  strains,  we  can  derive  the  following  expressions 

,  1  A  p 

€<■  =  §k=  -«,l  =  «!  +  (5) 

3  P 

and 

e«i  =  -  Ap/p  =  «-  +  2e,.  (6) 

If  A  pip  =  0  then  e,.  =  i.e.  the  axial  strain  becomes 

equal  to  the  effective  strain. 

Measurements  of  A  p/p  in  creep  experiments  show 
that  it  usually  lies  in  the  range  of  10' 3  to  10" 3  [1], 
Thus  in  ductile  materials  the  correction  from  Ap/p  in 
the  estimate  of  ee  in  equation  (5).  is  likely  to  be  small. 
In  long  term  creep-fracture  experiments  with  struc¬ 
tural  materials,  however,  the  apparent  tensile  strain  at 
fracture  can  be  small  as  1",,.  In  that  instance  the  cor¬ 
rection  becomes  large  and  it  should  be  implemented. 

Since  the  cavitation  strain  measures  the  total 
volume  of  all  the  cavities,  it  can  be  estimated  theoreti¬ 
cal!  v  The  task  becomes  particularly  simple  if  the  fol¬ 
lowing  assumptions  are  made:  (a|  that  the  shape  of 
the  grain  boundary  cavities  is  lenticular,  and  (b)  that 
the  number  of  cavities  remains  fixed.  The  first 
assumption  omits  the  fact  that  cavities  can  propagate 
like  cracks  under  certain  conditions  f  1 2],  and  the 
second  that  cavities  usually  increase  in  number  as 
well  as  in  volume  with  time  [  13)  Both  factors  can  be 
taken  into  account  by  increasing  the  complexity  of 
the  calculation 

We  assume,  as  illustrated  in  Fig  I.  that  cavities 
form  onlv  those  boundaries  which  are  approximated 
normal  to  rite  applied  tensile  stress,  that  their  average 
spacing  in  a  grain  facet  is  z,  and  that  their  shafic  is 
described  by  two  parameters  r  and  t)  From  this  geo¬ 
metry  a  relationship  between  and  •).  the  area  frac¬ 
tion  of  cavities  in  a  gram  facet  can  be  found  using  the 
expressions  in  Refs  [b.  '] 


Fh  =  n  sin3  0, 


(9| 


and  where  'J  is  the  grain  size. 

In  the  constrained  cavity  growth  mechanism,  the 
rate  of  growth  of  cavities  becomes  controlled  by 
matrix  creep.  Thus,  even  though  cavities  on  a  grain 
facet,  such  as  AB  in  Fig.  I,  grow  by  diffusion,  they  do 
so  under  a  tensile  stress  which  is  much  smaller  than 
the  applied  tensile  stress.  It  we  assume  this  local  stress 
to  be  <t'.  then  the  following  equation  can  be  written 
for  the  rate  of  growth  of  damage  A 


H)} 


where  c  is  a  normalized  unit  of  lime,  defined  as 
4  /I  32  F3  Vn  SDA 


where  t  is  real  time,  ft  is  the  atomic  volume  and  SDb  is 
the  boundary  width  times  the  grain  boundary 
diffusivity. 

Equation  (10)  is  a  constitutive  equation  for  de¬ 
scribing  the  rate  of  accumulation  of  damage  with 
time.  Interestingly,  it  is  universal.  The  influence  of  the 
microstructural  parameters,  such  as  A,  and  the 
material  parameters,  eg.  SDh.  and  the  stress,  have 
been  conveniently  normalized  with  respect  to  time  in 
£  defined  in  equation  (11).  The  plot  of  A  vs  £  is  given 
in  Fig.  2.  If  we  represent  this  master  curve  as  function 
of  £  such  that 

-4  =  fj(£)  (12) 

and  substitute  for  A  in  equation  (7),  we  obtain  the 
following  equation 


0-164 


Fig  3  A  plot  ol  t  s x  _  obtained  from  integrating  equation 
l  lill  Here  3  is  the  cavitation  damage  and  ^  is  a  normal¬ 
ized  unit  of  lime 


where 
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Fig.  3.  A  master  plot  showing  the  change  in  e„  with  =.  The 
microstructural  parameters,  loading  parameters,  and 
material  parameters  have  been  isolated  into  parameter  C 
[equation  (14)]  and  into  </r  [equation  (II)]. 


l°gioe«  —  C  +  (3/2)  log  10  (fif(c)!  (13) 


where 


„  ,  i  n  2  -  3  cos  6  +  cos 

C  =  lOg.o  <TZ 


12 


+  log.o 


(14) 


Thus  equation  (13)  can  be  used  to  compare  analysis 
with  the  measurement  of  the  change  in  density.  with 
time.  If  plotted  on  a  log-log  scale,  all  data  should 
overlap  on  a  single  master  curve  if  the  shift  factors  C 
and  c/f  are  adjusted.  Both  factors  depend  on  material, 
micrqstructurai,  and  loading  parameters.  C  depends 
on  0,  /.,  and,  as  given  in  equation  (14)  and  i/t  on  0, 

< x,  and  6Dh  as  described  by  equation  (II).  The  master 
curve,  which  is  exactly  defined  when  C  =0,  is  shown 
in  Fig.  3. 

The  validity  of  the  universal  curve  in  Fig.  3  is  nicely 
demonstrated  from  the  data  of  Needham  and  Green¬ 
wood  [1],  They  have  obtained  plots  for  log  e„  vs  lot  t 
for  various  levels  of  maximum  principal  stress  in  ten¬ 
sile  tests  done  with  a  superimposed  hydrostatic  press¬ 
ure.  When  the  data  in  Ref.  [I]  are  plotted  with  appro¬ 
priate  shift  factors  then  excellent  agreement  with  the 
master  curve  is  obtained  as  shown  in  Fig.  4. 

The  analysis  by  Harris  et  al.  [8]  differs  from  the 
one  given  here  in  two  respects-  First,  their  analysis  is 
based  upon  the  Hull  and  Rimmcr[l<).  II]  derivation 
in  which  the  applied  tensile  stress  acts  directly  on  a 
transverse  grain  boundary,  whereas  in  the  constrained 


cavity  growth  model  [3. 4]  the  local  stress  on  a  grain 
facet  may  be  lower  than  the  applied  stress,  its  magni¬ 
tude  depending  on  the  micromechanics  of  stress  redis¬ 
tribution  on  the  scale  of  the  grain-size.  In  the  present 
paper  this  effect  is  accounted  for  by  using  a  normal¬ 
ized  function  of  time.  The  second  difference  is  that  a 
master  curve  has  been  derived  for  the  change  in  rela¬ 
tive  density  with  time  (Fig.  3).  Convincing  evidence  is 
presented  for  the  validity  of  the  master  curve.  Several 
other  experimental  studies  and  a  phenomonological 
formulation  of  life  prediction  using  f„  as  the  critical 
parameter  [14]  have  been  reported  in  the  literature.  A 
review  of  these  results  is  not  attempted  here. 

Before  closing  this  section,  one  feature  of  equation 
(7)  is  worth  emphasizing  It  explains  how  the  change 
in  density  in  the  specimen  at  the  time  of  fracture  is 
related  to  microstructural  parameters  I).  '/.  and  y.  The 
relationship  between  for  different  values  of  0  is 
plotted  in  Fig.  5.  Note  that  for  a  given  damage,  say 
A  =  0.5.  the  change  in  density  increases  as  0  in¬ 
creases.  Also,  the  change  in  density  is  proportional  to 

✓./a 


EQUATIONS  FOR  FLOW 

The  next  step  is  to  consider  the  constitutive  equa¬ 
tions  for  flow  in  a  cavitating  solid.  In  particular  we 
wish  to  explain  how  cavity  growth  can  enhance  the 
deviatoric  component  of  the  strain-rate.  e,.  Dyson  et 
al.  [2]  have  suggested  that  the  constrained  cavity 
growth  mechanism  can  explain  this  effect.  Although 
similar  in  concept,  our  analysis  differs  from  theirs  in 
that  the  deformation  is  separated  into  two  scalar  par¬ 
ameters.  e„.  which  is  merely  the  fractional  change  in 
volume  ande,.  which  accounts  for  the  shape  change  at 
constant  volume.  Since  the  measured  principal  strain- 
rates  contain  a  component  of  volume  change,  the 
Levy-Mises  constitutive  equations  of  plasticity  [15], 
which  apply  only  when  deformation  occurs  at  con¬ 
stant  volume,  may  not  be  used  [2.  equation  (2)].  The 
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Fig  4  The  master  curve  m  Fig  5  shows  excellent  tit  with 
the  data  obtained  by  Needham  and  Greenwood  |  I  ] 
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Fig.  5.  A  plot  of  e„  vs  A.  as  prescribed  by  equation  (7). 

present  analysis  is  restricted  to  the  description  of  ie  in 
terms  of  the  principal  stresses. 

The  formulation  of  creep  enhancement  is  done 
through  a  new  description  of  the  deviatoric  stress. 
Normal  power  law  creep  is  represented  by  the  follow¬ 
ing  equation 

=  Vtc"  (15) 

where  ae ,  the  von-Mises  deviatoric  stress  has  the  form 

=  yiUfft  -  ff2)2  +  ( o2  -  <r3)2  +  (<r3  -  er,  )2 } . 

(16) 

In  extending  equation  (15)  to  the  situation  of  cavi¬ 
tation  enhanced  creep-rate,  some  features  of  the 
physical  mechanism  which  is  responsible  for  the  en¬ 
hancement  must  be  accounted  for.  For  example  in 
Fig.  1,  the  strain-rate  in  direction  1  caused  by  a  com¬ 
pressive  principal  stress  tr2  remains  unchanged 
because  cavity  growth  in  interface  AB  depends  only 
on  av  Also  for  the  same  reason  a  compressive  tr,  will 
not  produce  an  enhancement  in  the  rate  of  creep. 
Creep  is  enhanced  through  cavitation  only  when  the 
principal  stress  components  are  tensile.  Recognizing 
the  need  for  distinguishing  between  a  tensile  and  a 
compressive  principal  stress  wc  define  a  vector  6  such 
that 

<5,  =  +  1  if  a,  >  0 
=  -1  if  <t,  <0  (17) 

where  <Xj  are  the  principal  stresses  defined  by  the 
stress  tensor 


■»,  0  0  ' 

[a]  =0  <r,  0 

-0  0 

We  still  define  the  effective  stress  by  the  von-Miscs 
expression  but  replace  a,,  etc.  in  that  expression  by 
[1  +  ]x(l  +  <$,)](T,.  etc.  Thus  if  <5 ,  —  -  I  then  we 
recover  the  simple  von-Mises  form  but  if  <5,  =  +  1. 
then  a  stress  concentration  factor  of  (1  +  x)  is  ob¬ 
tained.  Substituting  these  into  the  von-Mises  equation 
gives  an  expression  for  the  deviatoric  stress.  which 
may  be  used  for  describing  cavitation  enhanced  creep 

qe  =v  Kid  +  i*)(<ti  -  Oi)  +  -  (r2f>2);2 

+  {(I  +  jx)(<7,  -  er3)  +  ]x(o202  -  rr3<f3 )]2 
+  !(1  +  33t)(fT3  -  <T, )  +  jx(rr3d3  -  a,S,  (j2]1  2. 

(18) 

Note  that  setting  x  =  0,  recovers  the  pure  von-Mises 
form. 

A  few  specific  applications  of  equations  (18)  are 
interesting.  First  we  consider  the  example  of  a  tensile 
test  superimposed  with  a  hydrostatic  pressure,  p.  Let 
<r.  be  the  tensile  stress  applied  to  the  specimen  inside 
the  pressure  vessel.  If  [a.  -  p)  >  0.  then  <r,  =  a.  —  p, 
a2  =  cr3  =  -  p  and  S2  =  I  and  S2  =  d3  =  - 1.  If. 
however,  (tr.  —  p)  >  0.  then  S ,  =  -I.  Substituting 
these  in  equation  (18)  we  obtain 

q€  =  a.  +  x(<7.  -  p)  if  ff5  -  p  >  0 

and 

qe  =  n.  if  a.  -  p  <0.  (19) 

In  a  torsion  test  <r,  =  a,  and  ct2  =  -a,  while 
er3  =  0.  Therefore  =  +  1,  S2  =  —  1  and  we  obtain 
that 

+  i“(*  +  2))1'2  %  aJV/3(l  +  $a)  (20) 

where  the  approximate  result  neglects  second  and 
higher  order  terms  in  x.  Equations  (19)  and  (20)  pro¬ 
vide  a  means  of  comparing  the  effective  stress  in  a 
tensile  and  a  torsion  test.  If  p  =  0.  then  the  stress 
concentration  factor  in  a  tensile  test  is  (1  +  x).  In  a 
torsion  test  is  (1  +  l/3x). 

The  cavitation  enhanced  effective  strain-rate  is  cal¬ 
culated  by  substituting  qe  instead  of  a ,  in  equa¬ 
tion  [15]. 

For  the  purpose  of  comparing  theory  and  experi¬ 
ment.  the  parameter  x  can  be  regarded  as  an  adjust¬ 
able  parameter.  The  equivalent  parameter  in  Dyson  et 
al.'s  analysis  is  A,:  and  they  are  related  through  the 
equation  (1  +  x)  =  1/(1  -  A,).  Upon  applying  their 
own  analysis  to  their  tension  torsion  experiments. 
Dyson  cl  «/.  found  that  -1 .  =  0.33  gave  satisfactory 
correlation  with  data.  Applying  it  to  the  superim¬ 
posed  hydrostatic  pressure  creep  experiments  by 
Needham  and  Greenwood,  a  value  A,  =  0.5  was 
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needed  to  give  good  fit  with  the  data.  The  application 
of  the  present  analysis  to  the  same  sets  of  data  shows 
that  using  x  =  0.42  (or  A ,  =  0.3)  gives  good  corre¬ 
lation  with  the  tension/torsion  data  and  z  =  1.5  (or 
A ,  =  0.6)  gives  satisfactory  correlation  with  Needham 
and  Greenwood's  data.  The  difference,  although  nu¬ 
merically  small  is  a  result  of  the  different  approach 
used  in  calculating  the  effective  strain-rate  from  the 
effective  stresst. 

An  extension  of  equation  (18)  when  the  loading  is 
cyclic  is  given  in  the  Appendix. 

MONKMAN-GRANT  DUCTILITY 

In  a  simple  uniaxal  test,  when  p  =  0.  the  minimum 
creep  rate,  and  time-to-fracture,  </.  are  often 
related  by  ictf  =  fmr  a  constant  in  long  term  fracture 
tests.  The  constrained  cavity  growth  model  success¬ 
fully  explains  this  result  [3. 4], 

If  the  tensile  test  is  done  in  a  superimposed  hydro¬ 
static  pressure  then  the  time-to-fracture  should  corre¬ 
late  with  the  maximum-principal  stress,  a,  since,  in 
the  constrained  cavity  growth  model,  the  fundamental 
mechanism  of  cavity  growth  is  diffusion.  The  Monk- 
man  Grant  equation  should  then  take  the  form 

kpa"ltJ=tm,  (21) 

In  tests  with  superimposed  hydrostatic  pressure 
<t,  =  <Tj  -  p.  Needham  and  Greenwood  did  in  fact 
find  that  (<r.  -  p)  and  tf  are  related  through  a  power- 
law  equation.  It  would  also  be  interesting  to  consider 
if  the  growth  in  cavitation  strain.  e„,  with  time  also 
follows  a  similar  rule.  A  logarithmic  plot  ofe„  normal¬ 
ized  with  respect  to  ( bz  -  p)*1,  against  time  is  shown 
in  Fig.  6.  The  results  correlate  quite  well  but  they 
separate  into  two  groups,  one  for  p  =  0.1  MPa  and 
3.5  MPa.  and  the  other  for  p  =  6.5  MPa,  10.4  MPa 
and  13.8  MPa.  A  possible  explanation  is  that  the 
cavity  spacing  for  these  two  groups  of  results  was 
different.  Possibly,  the  material  contained  two  types 
of  nucleation  sites,  one  requiring  a  greater  threshold 
stress  than  the  other.  Presumably  at  the  high  tensile 
stresses  (p  -  0.1  and  3.5  MPa)  both  types  nucleated 
leading  to  a  smaller  cavity  spacing  than  at  the  lower 
stresses  (p  >  6.5  MPa)  when  only  one  type  nucleated. 

Whereas  the  form  of  equation  (21)  has  a  physical 
basis  when  tests  are  done  in  uniaxal  loading  under  a 
superimposed  hydrostatic  pressure,  caution  must  be 
exercized  if  it  is  to  be  extended  to  triaxial  tensile  load¬ 
ing  of  an  isotropic  material.  The  Monkman  Grant 
equation  arises  from  the  fact  that  cavities  form  inho- 
mogeneously  in  the  material,  and.  therefore,  their 
growth  becomes  limited  by  compatibility  constraints. 
If  the  material  properties  were  ideally  isotropic,  i.e. 

t  For  example,  for  the  case  of  a  torsion  test,  equation 
(201  gives  the  cavitation  enhanced  effective  stress.  which 
replaces  a,  in  equation  (15)  to  give  the  cavitation  enhanced 
effective  strain-rate.  Dyson  cl  «/.  get  a  different  result  for 
the  effective  strain-rate  because  they  calculate  it  through 
the  Levy  Mises  equations  (their  equations  2  and  A 10). 
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Fig.  6.  A  normalized  plot  of  data  from  Ref.  [I]  showing 
the  correlation  between  d«„  dr  and  the  maximum  principal 

stress. 


the  grain-size  were  uniform  and  equiaxed  with  all 
boundaries  having  the  similar  density  and  types  of 
nucleation  sites,  then  under  equi-triaxal  loading  the 
compatibility  constraint  should  become  unimportant, 
and  the  simple  Hull  and  Rimmer  analysis [7.  10.  11.] 
should  apply.  In  the  author's  view  [4]  any  realistic 
situation  is  likely  to  deparl  considerably  from  the 
ideal  conditions,  and  equation  (2!)  should  continue  to 
hold  even  in  a  general  multia.xial  stress  state. 

DISCUSSION 

Three  inter-related  topics  have  been  discussed  in 
this  paper.  First,  the  need  for  separating  the  cavita¬ 
tion  dilation  strain  and  the  creep  strain  is  emphasized. 
Second,  it  is  shown  that  the  growth  of  cavitation 
strain  with  time  can  be  represented  on  a  master  curve 
provided  that  the  strain  and  the  time  are  normalized 
with  respect  to  the  microstructural  and  the  material 
parameters  such  as  0.  Q  and  <>£>,,.  The  master  curve 
can  be  used  to  predict  the  residua!  life  of  a  specimen 
after  it  has  suffered  some  cavitation  damage,  if  the 
normalization  parameters  C.  in  equations  (13)  and 
(14),  and  £/f,  in  equation  (II).  can  be  calculated.  The 
parameters  can  be  determined  experimentally  from 
the  knowledge  of  e/.  i.e.  the  volumetric  strain  at  frac¬ 
ture.  and  the  time-to-fracture.  it.  Substituting  .4  =  1 
in  equation  (13)  will  give  C  =  log10e^.  while  the  ratio 
i/t  may  be  determined  from  the  knowledge  that 
If  =  0.01 1  (see  Fig.  2)  and  by  using  c  t  =  if;tf.  Once 
the  shift  factors  are  established  the  master  curve  can 
be  used  to  predict  the  evolution  of  damage  with  time. 

The  other  topic  in  the  paper  deals  with  cavitation 
enhanced  deviatoric  creep.  The  peculiarities  of  the 
underlying  mechanism  for  this  phenomenon  are 
rather  novel  in  that  grain  facets  which  contain  grow¬ 
ing  cavities  behave  as  if  they  do  not  support  much 
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normal  traction,  even  though  the  grain  boundary  has 
not  physically  separated.  As  a  result  the  strain-rate,  in 
a  uniaxial  test,  is  enhanced  by  a  factor  U  +  xf  where 
a  is  a  factor  of  order  one.  and  where  n  is  the  power- 
law  stress  exponent.  The  strain-rate,  therefore,  may  be 
enhanced  by  almost  a  factor  of  10.  The  parameter  a 
has  been  left  as  an  adjustable  parameter  in  the  present 
analysis.  Physically  it  is  a  measure  of  the  extent  to 
which  the  normal  stress  is  transferred  from  the  grain 
facets  on  which  cavities  are  growing  to  the  adjacent 
surrounding  regions. 

Factors  which  enhance  the  strain  rate  arising  solely 
from  the  growth  of  cavities,  relative  to  the  rate  of 
creep  in  the  grain  matrix,  are  likely  to  increase  the 
value  of  x.  Since  cavities  must  grow  by  diffusion  to 
produce  this  effect,  while  the  grain  matrix  is  assumed 
to  deform  by  power-law  creep,  the  following  factors 
are  likely  to  increase  the  value  of  x:  smaller  '/.  because 
that  leads  to  faster  cavity  growth  by  reducing  the  dif¬ 
fusion  distance:  smaller  applied  stress  because  dif- 
fusional  cavity  growth  varies  linearly  with  stress  while 
creep  varies  as  stress  to  a  power  of  4  to  5;  and  inter¬ 
mediate  temperatures  because  diffusional  growth  is 
usually  dominated  by  grain  boundary  diffusion  while 
power-law  creep  is  controlled  by  lattice  diffusion.  The 
enhancement  in  the  creep-rate  would  also  depend 
upon  the  spatial  redistribution  of  the  load  which  is 
shed  by  the  grain  facets  and  transferred  to  the  sur¬ 
rounding  grain  matrix.  In  a  material  with  high  strain 
rate  sensitivity  the  load  will  quickly  spread  across  the 
load  bearing  regions  and,  therefore,  lead  to  a  greater 
increase  in  the  measured  creep-rate. 

Another  interesting  point  about  the  above  mechan¬ 
ism  is  that  the  creep-rate  is  insensitive  to  the  size  of 
the  cavities,  the  only  necessary  condition  is  that  cavi¬ 
ties  should  have  nucleated  and  should  be  growing. 
This  result  was  confirmed  by  Needham  and  Green¬ 
wood  who  found  the  creep-rate  did  not  depend  on  the 
extent  of  prior  damage,  and  by  Dyson  et  al.  who 
found  that  the  effect  occurred  only  if  the  cavities  have 
been  pre-nucleated.  This  raises  an  interesting  appli¬ 
cation  of  the  phenomenon.  Often  in  engineering 
materials,  there  is  a  need  to  determine  when  the  cavi¬ 
ties  nucleate  in  a  stress-rupture  test.  If  an  alternate 
compressive  and  tensile  load  were  applied  to  the 
specimen,  then  the  compressive  creep-rate  and  the 
tensile  creep  rate  should  remain  equal  until  the  cavi¬ 
ties  nucleate.  Thereafter,  the  tensile  creep-rate  will 
become  significantly  faster  than  the  compressive 
creep-rate. 

SUMMARY 

1.  In  long  term  experiments  with  structural 
materials,  when  the  fracture  strain  is  small,  the  strain 
due  to  cavitation  (change  in  volume)  and  the  strain 
due  to  creep  (shape  change  only)  must  be  separated. 
The  cavitation  strain.  c0.  defined  in  equations  (2)  and 
(4)  and  the  creep  strain,  e,.  in  equation  (3),  are  inde¬ 
pendent  quantities  in  any  strain-state. 


2.  The  constitutive  flow  equation  for  e,„  however, 
may  be  coupled  to  and  enhanced  by  cavitation  if  the 
cavities  are  growing  by  the  constrained  cavity  growth 
mechanism.  Then  the  effective  stress  and  «,  are  still 
related  through  the  usual  (low  equation,  but  the  effec¬ 
tive  stress.  oe,  is  no-longer  described  by  the  customary 
von-Mises  equation. 

3.  A  possible  new  definition  for  the  deviatoric 
stress.  qe.  is  given  in  equation  (18).  It  depends  not  only 
on  the  principal  stresses  but  also  on  the  sign  of  the 
principal  stresses,  since  a  tensile  stress  is  required  to 
grow  cavities  by  diffusion.  The  expression  contains  an 
adjustable  parameter  x. 

4.  x  can  be  determined  by  comparing  the  flow  be¬ 
havior  in  simple  tension  and  simple  compression. 

5.  The  relative  rates  of  flow  in  tension  and  com¬ 
pression  may  also  be  used  to  evaluate  cavity  nu- 
cleation,  since  creep  enhancement  in  tension  is  poss 
ible  only  after  the  cavities  have  nucleated. 
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APPENDIX 

In  a  fatigue  experiment  where  compressive  deformation  is 
imposed  after  allowing  the  cavities  lo  grow  in  tension,  some 
enhancement  in  the  compressive  creep-rate  may  also  be 
expected  because  the  cavities  sinter  in  compression  which 
will  lead  to  a  displacement  rate  of  A  and  B  toward  each 
other  (sec  Fig.  1).  The  enhancement  factor  in  compression 
is  likely  to  be  smaller  than  in  tension  Let  the  corresponding 
value  of  x  in  compression  he  x'.  then  the  principal  stress 
in  the  von-Mises  equation  should  be  written  as 
[1  +  }x(l  +(),)+  jx'd  -  4,)]<r,.  etc  Note  that  they 
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reduce  to  [1  +  $*(1  +  <5 ,  >]  .  and  so  on.  for  tension  and  +  {(1  +  +  Jx'Kfjj  ~  fli> 

[I  +  Ja'(l  -  <),)]<t,  for  compression.  With  these  substitu-  . 

tions  the  expression  for  von-Mises  effective  stress  becomes:  +  -  x  -  ''s'Gll 

=v  KlO  +  1*  +  iar'Ker,  -  <r,|  +  !(l  +  +  }x')(o3  -  <T|  I 


+  -  or" )(<>, <r,  -  d2<T2)|J 


+  j(i  -  x')(dji73  -  d|(7.  If-2]'  J.  (Al) 
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Abstract  Gaseous  species  can  diffuse  along  i  boundary  and  oxidize  impurities  or  precipitates  In  nickel 
alloys  O,  can  attack  C  or  carbides  to  produce  CO  at  high  virtual  pressures.  This  provides  a  driving 
force  for  cavity  nucleation  and  growth  at  particle  interfaces,  which  acts  in  addition  to  the  applied  stress 
in  forcing  fracture  In  this  paper  it  is  shown  that  the  oxidation  pressure  can  become  comparable  to.  or 
even  larger  than,  the  stresses  normally  applied  in  stress-rupture  tests.  The  dependence  of  this  pressure  on 
temperature  and  oxygen  pressure  is  derived  for  various  oxidation  reactions.  Since  oxide-dispersion- 
strengthened  alloys  arc  known  to  be  more  susceptible  to  cavitation  than  precipitate  strengthened  alloys, 
it  is  also  possible  that  internal  oxidation  enhances  cavitation  by  producing  oxide  particles  from  elemen¬ 
tal  impurities,  or  by  converting  metal  carbide  and  metaUintermetallic  interfaces  to  metal-oxide  inter¬ 
faces.  The  kinetic  steps  which  may  become  rate  limiting  when  cavity  nucleation  and  growth  is  induced 
by  oxidation,  are  discussed 


Resume  Les  gaz  peuvent  diffuser  le  long  d'un  joint  de  grains  et  oxyder  des  impuretes  ou  des  precipites. 
Dans  les  alliages  de  nickel.  1'  oxygene  02  peut  attaquer  le  carbone  ou  les  carbures  et  former  de  I'oxyde 
de  carbone  CO  sousde  fortes  pressions  virtuelles.  Ceci  induit  une  force  motrice  pour  la  germination  et  la 
croissance  des  cavites  aux  interfaces  avec  les  particules.  force  mortice  qui  s'ajoute  a  la  contrainte 
appliquee  pour  provoquer  la  rupture.  Dans  cet  article,  nous  montrons  que  la  pression  d'oxydation  peut 
devenir  comparable  ou  mime  superieure  aux  contramtes  appliquees  normalement  au  cours  d'essais  de 
rupture  sous  contrainte.  Nous  avons  obtenu  la  variation  de  cette  pression  en  fonction  de  la  temperature 
et  de  la  pression  d'oxygene  pour  diverses  reactions  d'oxydation.  Les  alliages  renforccs  par  une  dispersion 
d'  oxyde  presentant  plus  frequemment  le  phenomene  de  cavitation  que  les  alliages  renforces  par  des 
precipites,  il  est  egalemcnl  possible  que  I'oxydalion  interne  favorise  la  cavitation  cn  transformant  des 
impuretes  sous  forme  d'elements  en  particules  d'oxyde.  ou  en  transformant  les  interfaces  metal-carbure 
et  metal  compose  intermetallique.  cn  interfaces  metal  oxyde  Nous  disculons  les  elapes  de  la  cinctique 
susceptibles  de  limiter  la  vitesse  lorsquc  la  germination  et  la  croissance  des  cavites  provienent  de 
I'oxydalion. 

Zusammenfassung— Gasfbrmige  Stoffe  konnen  entlang  Korngrenze  diffundieren  und  Verunreinigungen 
Oder  Ausscheidungen  oxidieren.  In  Nickel-Legierungen  kann  02  bei  hohem  virtuellem  Druck  mit  C 
oder  Karbiden  zu  CO  reagieren.  Dadurch  konnen  Hohlraume  gebildet  werden  und  an  den  Ausschei- 
dungsgrenzfliichen  wachsen.  Dieser  Vorgang  unterstiitzt  die  angelegte  Spannung  bei  der  Bruchbildung. 
In  dieser  Arbeit  wird  gezeigt,  daB  der  Oxidationsdruck  vergleichbar  oder  sogar  groOer  werden  kann  aid 
die  Spannung,  die  iiblicherweise  in  Spannungsbruchversuchen  angewendet  wird.  Die  Temperatur-  und 
Sauerstoffdruckabhangigkeit  dieses  Druckes  wird  fur  verschiedenc  Oxidationsreaktionen  abgeleitet.  Es 
ist  bekannt,  daB  Legierungen,  die  mit  Oxidpartikeln  gehartet  sind.  starker  zur  Hohlraumbildung  neigen 
als  solche,  die  mit  Ausscheidungen  gehartet  sind.  Daher  ist  es  auch  moglich,  daB  innere  Oxidation  die 
Hohlraumbildung  dadurch  fordert,  daB  elementar  vorliegende  Verunreinigungen  zu  Oxidteilchen  oder 
Metailkarbiden  und  Metall-intermetallische  Grenzflachen  zu  Melall-Oxidgrenzflachen  umgewandelt 
werden.  Die  kinetischen  Schrilte,  die  bei  Hohlraumnukleation  und  -wachstum  durch  Oxidation  gesch- 
windigkeitsbestimmend  werden  konnen.  werden  behandelt. 


1.  INTRODUCTION 


Nearly  always,  engineering  materials  face  aggressive 
environments  during  service,  which  is  why  a  funda¬ 
mental  understanding  of  phenomena  such  as  hydro¬ 
gen  and  liquid  metal  embrittlement  and  stress-corro¬ 
sion  cracking  is  so  important.  At  high  temperature, 
the  sensitivity  of  Ni  alloys  to  air-environment  has 
been  reported  by  Sadananda  et  at.  f  1.  2.  3],  Recent  ex¬ 
periments  by  Woodford  and  Bricknell  [4,  5. 6]  have 


established  beyond  doubt  that  oxygen  is  an  aggressive 
species  in  the  creep-fracture  behavior  of  nickel  base 
alloys.  They  have  shown  [5]  that  cavities  can  be 
induced  in  Ni  even  when  it  contains  small  amounts  of 
C  ( —  100  ppm),  and  at  oxygen  pressures  which  are  low 
enough  to  be  reducing  with  respect  to  NiO.  They 
propose  that  the  cavities  arc  bubbles  of  CO.  This  re¬ 
inforces  this  author's  belief  that  oxygen  attack  in 
alloys  of  Ni  does  not  involve  a  new  mechanism  of 
creep  fracture,  but  that  it  accentuates  the  usual  creep- 
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fracture  mechanism,  i.e.  cavity  nucleation  and  growth, 
which  occurs  even  in  an  inert  environment  when  a 
tensile  stress  is  applied  This  intuition  is  based  on  the 
study  of  phenomenogical  data  which  compares  frac¬ 
ture  in  vacuum'  and  air  environments,  eg  Refp]  In 
general,  one  observers  (a I  that  althouvt'  fiaclure 
occurs  sooner  in  air  than  in  vacuum,  cavitation  i*  still 
the  principal  feature  of  fracture  surfaces,  and  ibl  hat 
strong  nickel  base  alloys  are  often  more  notch  sensi¬ 
tive  in  air  than  in  vacuum  although,  again,  cavitation 
is  the  cause  of  fracture.  The  second  point  car,  be 
rationalized  as  being  due  to  enhanced  cavi.jiion  by 
the  diffusion  of  oxygen  from  the  surface 

In  this  paper,  therefore,  the  earlier  work  on  nu¬ 
cleation  and  growth  of  cavities  under  stress  [8.  9.  10], 
is  extended  to  include  the  additional  thermodynamic 
driving  force  for  cavitation  due  to  the  gas  produced 
by  the  reaction  between  O  and  C  (and  its  com¬ 
pounds).  Like  the  tensile  stress  across  the  boundary, 
the  gas  pressure  attempts  to  pry  the  boundary  apart. 
The  kinetics  of  such  a  process,  however,  is  much  more 
complex,  and  can  be  considered  only  qualitatively. 
For  example,  trace  impurities  can  influence  the  dif¬ 
fusion  rate  of  O  through  the  interface,  and  surface 
structure  can  influence  the  dissociation  reaction  of 
molecular  oxygen.  One  must  also  think  about  the 
mobility  of  CO  molecules  in  the  interfaces,  since  a 
boundary  which  is  intact,  these  must  diffuse  and  co¬ 
alesce  or  bind  with  vacancies  to  produce  a  bubble 
which  is  larger  than  a  critical  size. 

The  analysis  developed  here  is  general  and  may  be 
applied  to  other  reactions  such  as  sulphides  reacting 
with  02  to  produce  S02.  or  with  H2  to  produce  H2S, 
silicides  reacting  with  02  to  produce  SiO,  and  H2 
with  carbides  to  give  CH*.  Since  we  are  concerned 
here  with  nickel  alloys,  primarily,  we  shall  consider 
only  the  following  reactions 

Cr23C#(s)  +  302(</)  =  23Cr(s)  +  6CO (<?)  (Rl) 

2Cr7C3(s)  +  30 2(g)  =  14Cr(s)  +  6CO{g)  (R2) 
Cr3C2(.s)  +  O  2(g)  =  3Cr(s)  +  2CO(g )  (R3) 

2CO(g)  +  02(g)  =  2C02Uj)  (R4) 


% 


y, 

■i- 


where  the  letters  in  brackets  stand  for  solid  (s)  and  gas 
((/).  The  products  in  reactions  (Rl),  (R2)  and  (R3)  can 
react  further  with  02.  Carbon-monoxide  can  oxidize 
into  C02,  as  given  by  reaction  (R4).  In  fact,  at  tem¬ 
peratures  of  interest  and  at  p„,  =  1  (where  p  is  the 
equilibrium  partial  pressure),  oxidation  to  C02  is  a 
favored  reaction,  i.e.  pf0j  "ill  be  greater  than  p(0. 
although,  kinetically,  it  may  be  more  difficult  to 
obtain  C02  because  the  reactions  must  occur  in  the 
solid  state.  Similarly,  Cr  may  oxidize  into  Cr20, 
which  will  also  increase  the  CO  pressure.  In  nickel 
base  superalloys,  the  presence  of  excess  Cr  would  sug¬ 
gest  that  reactions  (Rl)  to  (R3)  rather  than  those  pro¬ 
ducing  Cr203,  are  applicable.  Finally,  for  those 
instances  when  free  C  is  present  we  consider  the  reac¬ 
tion: 


2C  -r  O,  -  2CO  (R5) 


In  the  Inflowing  sections,  thermodynamics  and  kin¬ 
etics  are  considered  separately  The  discussion  con¬ 
centrates  more  on  nucleation  rather  than  on  growth 
because  it  is  felt  that  the  production  rate  of  CO  or 
CO,,  which  may  be  limited  by  the  diffusional  trans¬ 
port  ol  O  from  the  surface,  may  not  be  fast  enough  to 
keep  up  with  the  growth  of  cavities  when  they  grow 
large,  but  may  have  significant  influence  when  cavities 
are  of  nm  dimensions.  The  treatment,  however,  fully 
includes  the  issues  related  to  both  nucleation  and 
growth 


! 

i 
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2.  THERMODYNAMICS 

2.1  Gu\  pressure  us  u  driving  Jorte  for  nucleation  and 
growth 

Consider  a  cavity  in  a  grain  boundary.  A  normal 
traction.  oT,  is  applied  across  the  boundary  and  the 
cavity  is  pressurized  with  a  gas  at  pressure,  p.  Assume 
now  that  the  cavity  grows  by  a  volume  61'  by  the 
transport  of  matter  from  the  cavity  surface  to.  and 
into,  the  grain  boundary.  In  such  case  the  elastic 
energy  may  be  neglected  [8]  and  the  change  in  energy 
of  the  system.  will  involve  the  work  done  by  <rx, 
by  p,  and  by  the  'sintering  pressure',  2y/r  where  y  is 
the  surface  energy  and  r  is  the  radius  of  curvature 
(assuming  that  both  radii  of  curvature  in  the  surface 
are  equal).  Thus 


where  5%  is  the  small  normal  surface  displacement, 
and  the  integral  represents  integration  over  the  sur¬ 
face  S.  If  p  P  (2y/r)  or  if  the  cavity  is  of  an  equilibrium 
shape  so  that  r  is  uniform  at  all  parts  of  the  cavity 
surface  then  (p-2y/>)  becomes  constant  and  (3)  reduces 
to 

sw  =  (p -?y  +  o.'jsv.  (7) 

It  follows  that  in  calculations  of  cavity  nucleation 
and  growth  trx  may  be  replaced  by  (cr*  +  p  -  2y />) 
provided  that  the  cavity  has  a  near  equilibrium  shape 
or  if  (it,  +  p)  >  2y/r.  Earlier,  this  author  has  pro¬ 
posed  a  diffusive  model  for  cavity  nucleation  [10] 
where  a  cavity  grows  by  clustering  of  vacancies  until 
the  cavity  becomes  super-critical.  For  sub-critical 
voids,  the  surface-energy  factor  is  dominant  and  the 
cavity  grows  slowly;  in  this  instance  equation  (7) 
should  apply  because  the  cavity  is  small  and  slow 
growing  so  that  it  would  retain  its  equilibrium  shape. 
When  the  cavity  is  much  larger  than  the  critical  size 
then  the  2 y/r  term  should  become  relatively  small  and 
again  equation  (7)  should  apply,  provided  that  Sy  is 
independent  of  S.  Caution  is  needed  when  the  cavity 
is  just  a  bit  bigger  than  its  critical  size. 
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In  equation  (7).  p  is  the  mechanical  pressure  of  the 
gas.  If  the  gas  is  produced  by  chemical  reactions  (R1 
thru  R 5).  then  the  virtual  pressure,  also  p.  will  be 
given  by: 

f=lp  |K| 

where  /  is  the  activity,  or  the  fugacity.  and  /  is  the 
activity  coefficient.  /.  thus,  is  determined  through  the 
free-energy  of  the  reaction  in  question,  and  becomes 
the  ‘chemical-  driving  force  for  cavity  nucleation.  We 
now'  discuss  the  magnitudes  for  /  and  p 

2.2  fco  Produced  b  t  relictions  between  C.  carbides  and 
02 

The  fugacity  of  CO  for  reactions  (Rl  to  R5)  is  cal¬ 
culated  in  this  section.  The  standard  free  energy  of 
formation  of  the  different  carbides  of  Cr.  CO  and 
C02,  and  Cr203.  are  listed  in  Table  1.  The  values  for 
the  carbides  have  been  taken  from  Ref.  [12]  because 
they  were  obtained  by  seeking  equilibrium  with  CO 
gas  and  because  they  are  in  agreement  with  other 
measurements  where  equilibrium  with  CO  was  also 
used  to  obtain  thermodynamic  data  [14].  These 
numbers  differ  slightly  from  those  given  in  Ref.  [13] 
(by  less  than  .V’,,  for  Cr23C„.  4°„  for  Cr-C3  and  8°,, 
Cr3C2).  For  a  recent  review  see  Ref.  [15]. 

The  fugacity  of  CO  is  plotted  as  a  function  of  tem¬ 
perature  and  p o,  for  reactions  of  O,  with  Cr23C„. 
Cr,Cj.  Cr3C2  and  C.  in  Fig.  1.  Note  that  Cr23C6 
gives  the  smallest  values  for  the  fugacity  and  C  the 
highest.  Typically,  at  1000  K.  the  values  for  Cr23Ch 
are  about  10'  at  />„,  =  1.  and  103  at  p„,  =  10' 8.  For 
pure  C  the  fugacitics  are  a  factor  of  3  x  10J  higher,  at 
1000  K. 

At  this  point  it  would  be  interesting  to  estimate  the 
virtual  CO  pressure  in  experiments  of  Bricknell  and 
Woodford  [5],  Their  technique  for  controlling  /\, 
was  to  pack  the  Ni-specimens  which  contained  some 
C  in  Ni:NiO  powder  at  1273  K,  which  fixed  the  02 
activity  at  p0j  =  3  x  10“  ",  and  then  see  if  bubbles 
were  produced  in  the  boundaries  and  they  were. 
Assuming  a  carbon  activity  of  one,  we  obtain  from 
Fig.  1,  that  Jen  =  7  x  103.  If  instead  vie  assume  that 
at,  the  carbon  activity,  was  equal  to  the  average  C 
concentration  (100  ppm  by  weight).  then 
ac  =  5  x  10'*,  which  reduces  /,,,  to  3.4.  At  these  low 
pressures  once  may  assume  the  F  =  1.  which  would 


mean  an  effective  CO  pressure,  p, „  =  0.34  MPa 
(using  equation  8).  Since  cavities  were  observed  under 
these  conditions,  one  may  infer  that  a,  >  5  x  10  4 
because  of  non-equilibrium  segregation  in  the  boun¬ 
daries.  There  is  also  the  possibility  that  CO;  instead 
of  CO  was  being  produced.  That  reaction  at  1273  K 
and  pu.  =  3  x  10  4,  yields  p,„,  =  21  MPa.  which  is 
a  stress  large  enough  to  cause  cavitation  in  pure 
nickel.  The  enhancement  factor,  when  CO:  rather 
than  CO  is  the  reaction  product,  is  given  in  Fig.  2, 
and  was  obtained  by  considering  reaction  |R4| 

The  assumptions  made  in  the  calculation  of  curves 
in  Figs  I  and  2  need  to  be  re-emphasized  The  activi¬ 
ties  of  Cr.  C  and  the  carbides  have  been  assumed  to 
be  unity.  While  this  is  probably  true  for  the  carbides, 
since  they  do  not  dissolve  significant  amount  of  im¬ 
purities.  the  activities  of  Cr  and  C  would  depend  on 
the  alloy  composition  and  in  whether  or  not  the 
boundary  segregation  of  these  impurities  is  in  equilib¬ 
rium  with  the  matrix.  Decreasing  C  activity  decreases 
/,(>.  while  decreasing  Cr  activity  increases  /,,,.  In  case 
of  Cr23Ch.  /<„  approximately  increases  as  l.  ti,  ,  [4], 
Possibly  a  more  accurate  estimate  of  activities  in  real 
alloys  can  be  done.  but.  in  the  author's  view  a  system¬ 
atic  experimental  study  on  the  ternary  Ni  Cr-C  sys¬ 
tem  is  first  necessary  before  attempting  to  apply  the 
analysis  to  more  complex  alloys.  Another  assumption, 
that  Cr  rather  than  Cr203  is  the  reaction  product,  is 
justified  on  the  basis  that  excess  Cr  is  usually  present 
in  the  nickel-base  alloys.  However,  for  sake  of  com¬ 
pleteness.  analytical  expressions  for  /,,,  for  this  and 
the  other  reactions,  are  given  in  Table  2. 

The  fugacity  of  gases.  /,  can  be  related  to  their  par¬ 
tial  pressure,  p.  through  the  expression  in  equation  (8|. 
For  low  values  of  f  up  to  a  few  atmospheres,  the 
behavior  should  remain  ideal,  and  the  activity  coeffi¬ 
cient  F  =  1.  In  the  following  section  we  describe  a 
method  for  obtaining  approximate  values  of  F  at 
large  partial  pressures  of  the  gas. 

3.  FUGACITY  OF  NON  POLAR 
GASES  AT  HIGH  TEMPERATURE 
AND  PRESSURE 

3.1.  Etpmtion  of  state 

At  high  pressures  the  compressibility,  and.  there¬ 
fore.  the  free  energy  of  gases  becomes  pressure  depen- 


Table  1.  Free  energy  of  formation 


Standard  free  energy 

of  formation  Temperature 

AG„  range 


Compound 

AG„ 

{J  mole  ') 

Reference 

K 

COlg) 

-III  6  x  10' 

-87.57T 

11 

298  2300 

CO;  Ig) 

-393.8  x  10' 

—  0.84T 

11 

298  2000 

Cr;,C„M 

-3219  x  10' 

-76  5T 

12 

1150  1 31X1 

Cr-C ,  |sl 

-  147  1  x  10' 

-36  37T 

12 

11(8)  1720 

Cr jC.  is) 

-68  6  x  10' 

-  I8  4T 

12 

13(8)  15(8) 

Cr;0,  |sl 

-  1 121  8  x  It)' 

+  252.8T 

13 

8(8)  1600 

VV 


TC0; 
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Fig.  2.  Relative  fugacity  of  CO,  and  CO  for  oxidation  of 
CO. 


give: 


In  / 


V 


7*f>" 

nRT 
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Maron  and  Turnbull  argue  that  equation  lltll  should 
be  universal  if  p  and  T  are  normalized  with  respect  to 
the  critical  values.  />.  and  Tt.  in  van-der  Waal's  equa¬ 
tion  such  that 

P  .  T 

pr  =  -  and  7,  =  1 1 1 1 

Pc  T 

where  pr  and  7,  are  the  reduced  or  normalized  quanti¬ 
ties.  With  these  substitutions  HOI  becomes 


where 


In  /  '  =  \  an 

n  l 


Pr 

nRT, 


a„ 


x„n,  I 


P " 

T, 


1 1  2l 


1 1  Tt 


Maron  and  Turnbull  have  shown  that  the  above 
parametric  description  fits  the  P.i.T  data  to  an  accu¬ 
racy  of  better  than  5"„.  The  coefficients  «„  are  listed  m 
Table  3. 


Table  2.  Equilibrium  values  of  CO  fugacity  for  different  oxidation  reactions 


CO(i/)  +  }02(</l  =  C02(g) 

Cr2jC6(.v)  +  3O2(0)  =  23Cr(s)  +  6CO(y)  fco  =  8.2  x  10 


~  =  Po7  x  2.9  X  10” 5  exp ( 282.2  kj  mole  '  1  RT  ) 
fco 


5  exp(58  kJ  mole  1  R 7  1 

b"2u'  3 

2Cr,C,(s)  +  302(a)  =  14Cr(s)  +  6CO(y)  /co  =  8.9  x  103  exp(62  6  kj  mole ” '  RT I 

“ c , 

Cr3C2(s)  +  02(g)  =  3Cr(s)  +  2CO(g)  fco  =  1.3  x  10*  exp(77.3  kJ  mole  1  R  7  | 


C(s)  +  JO  2{g)  =  CO(</) 


a 

fco  =  3.8  x  10*  pl,2ar  exp(  1 1 1  6  kj  mole  1  RT\ 


4Cr2JC6  +  8102  -  46Cr202  +  24CO  /co  =  3.9  x  10  22  —  T-'”-1*  exp(2208  kj  mole  ‘ 1  KTl 


Table  3  Numerical  expressions  for  the  coefficients  u„  in  equation  (131  These  numbers  are  mil  dimensionless,  they  have 
been  calculated  such  that  R  =  0082054  bier  bar  K  1  mole”  1  and  p  and  Fare  expressed  in  bars  and  liters  in  equations 

(9MI3) 


u  =  A„  +  V  Am  TT 


*0 

A , 

Ai 

A  a 

A. 

Ac 

•c 

1.0196  x  10  3 

-2.1420  x  10  1 

0 

-3.2548  x  10  ; 

0 

0 

0 

a. 

0 

0 

3.6991  x  10  3 

0 

-4  3022  x  10' 

0 

1  S.  M  ■ 

10  •' 

0 

0 

-  1.251  x  It)  * 

0 

4  6408  *  10  4 

0 

-  1  SS7V  . 

|n  ‘ 

<*A 

0 

0 

1  7019  x  10  * 

0 

-  1  1221  x  10  ' 

0 

3  2830  x 

10  ‘ 

1264 
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Fig  3.  Variation  of  /.  the  activity  coefficient,  with  p  for 
CH4  and  CO  The  CO  data  is  from  Refs  [19  and  20]  and 
the  CFU  data  from  Ref.  [IS] 


3.2  Applications  of  equation  (/ 3)  to  carbon-monoxide 
qas 

Compressibility  data  for  CO  extends  up  to  1200 
bars  and  673  K.  Since  the  fugacities  calculated  for  car¬ 
bide  oxidation  reactions  (Fig.  1)  can  be  up  to  x  I03 
higher,  we  must  consider  the  suitability  of  extrapo¬ 
lating  equation  (12)  to  higher  pressures.  Since  that 
polynomial  contains  terms  with  Tr  and  pr  raised  to  a 
high  power,  small  errors  in  the  numerical  coefficients 
will  lead  to  large  discrepancies  in  the  estimate  of  F. 
To  reduce  this  uncertainty  we  consider  only  the  first 
order  term  in  equation  (12).  Thus 

In  F=A,Pr  (14) 

T, 

where  Ar  is  a  universal  dimensionless  constant  which 
must  be  determined  from  data.  Keeping  in  mind  that 
the  critical  parameters  for  CO  and  CH4  are 

CO:  pc  =  35.0  bars  7]  =  134.18  K  ( 

CH4  p,  =  45.8  bars  Tc  =  190.68  K 


we  obtain  the  plots  show  n  in  Fig.  3  for  the  fugacity 
data,  according  to  equation  (14).  The  fit  with  (14)  is 
quite  good  when  pr  <  75,  but  above  that  the  value  of 
F  given  by  (14)  overestimates  the  data.  By  extrapola¬ 
tion.  we  may  assume  that  the  value  of  r.  if  calculated 
from  equation  (14).  will  always  be  >  to  the  actual 
value.  If  this  value  is  substituted  in  equation  (8),  then 
we  shall  obtain  the  lowest  possible  value  for  p. 

The  universality  of  equation  (14)  is  checked  by  cal¬ 
culating  A,  tor  CO  and  CH„  at  different  temperatures. 
These  results,  on  the  basis  of  the  straight  lines  drawn 
in  Fig.  3.  are  tabulated  in  Table  4.  A  mean  value  of 
Ar  =  8.57  x  10  2  is  obtained  with  a  maximum  devi¬ 
ation  of  less  than  7l’n. 

The  equilibrium  values  for  pt(>  are  obtained  by 
applying  equations  (14)  and  (8)  to  the  fugacity  curves 
in  Fig.  1.  The  results  for  the  reactions  (Rl)  and  (R5) 
are  shown  in  Fig.  4.  Remember  that  these  are  subject 
to  the  assumptions  discussed  in  this  section.  In  gen¬ 
eral  we  expect  that  the  numbers  are  quite  good  as 
long  as  /v<,  is  less  than  3  x  I03  bars.  When  the  esti¬ 
mated  /Vo  is  greater,  then  the  real  values  are  likely  to 
be  even  higher  than  those  derived  in  Fig.  4.  Neverthe¬ 
less  results  in  Fig.  4  provide  some  interesting  insights. 
When  pn,  =  1.  the  CO  pressure  will  be  greater  than 
103  MPa.  a  value  which  is  even  larger  than  the  stress 
applied  in  most  stress-rupture  tests.  The  pressure  falls 
sharply  when  the  pn,  drops  below  10  “  bars.  As  a 
guideline,  one  can  expect  that  when  pn,  is  below 
10  10  bars,  the  oxidation  reaction  would  be  thermo- 
d  >amics  limited,  while  above  10  “  bars  the  kinetics 
of  the  reaction  may  be  more  important  Also,  for  a 
fixed  value  of  /v,  the  environmental  attack  should  be 
most  severe  at  intermediate  temperatures;  at  high 
temperatures  pco  becomes  small,  while  at  the  lower 
temperatures  the  kinetics  becomes  slow  The  question 
of  kinetics  is  discussed  in  the  next  section 

4.  KINETICS 

4.1  The  heteroqeneous  mteleation  model 

The  mechanism  of  cavity  nucleation  in  creep-frac¬ 
ture  remains  controversial.  Limited  observations  sug¬ 
gest  that  cavities  do  not  nucleate  below  a  threshold 
stress  [21],  that  they  increase  in  number  with 
time  [22]  and  that  triple  junctions  formed  by  second 
phase  particles  and  grain  boundaries  are  the  preferred 
sites  for  nucleation.  There  is  more  in  the  literature 


Table  4  A  tabulation  of  the  constant  A,  in  equation  1141  based  upon 
thermodynamic  data  for  CO  and  CH4 


T 

K 

p  (range) 
bars 

Tr 

p,  (range) 

A. 

CO 

673 

25  - 1 2(H) 

5.02 

0.7-34.3 

9.18  x  10  ! 

CO 

573 

25  1200 

4.27 

0.7  34  3 

8,46  x  10  2 

CH4 

i.  3 

2000  85(X) 

3.53 

43.7  185.6 

8,54  x  10  ’ 

ch4 

573 

2000  85(H) 

301 

43.7-185.6 

8.34  x  10  1 

ch4 

473 

2000  85(H) 

2.48 

43.7  185.6 

8.35  x  10  : 

Average  value 

8.57  x  10 

RAJ  INTERGRANULAR  CREEP  FRACTURE  IN  AGGRESSIVE  ENVIRONMENTS 


1265 


'4.' .  «  « «  > 


I03/T  (K ) 

fig  4  Carbon-monoxide  partial  pressure  for  two  oxi¬ 
dation  reactions,  it,  i  carbon  activity. 


than  cited  here,  but  still  further  work  is  needed  before 
the  picture  becomes  clear.  For  the  present,  we  shall 
assume  that  the  model  published  earlier  [10]  is  viable 
The  cavities  can  nucleate  under  applied  stress  by  a 
mass  transport  process  provided  that  the  barrier  for 
cavity  nucleation  due  to  surface  energy,  is  not  very 
large  A  useful  physical  parameter  to  characterize  the 
nucleation  process  is  an  incubation  time  which  marks 
the  transition  fn  extremely  slow,  sub-crilica! 

growth  ol  cavities  to  fast,  super-critical  growth  We 
can  appreciate  that  the  incubation  time  will  depend 
upon  the  volume  of  the  c...  of  critical  size 

(expressed  in  unils  of  atomic  vo  i.  and  upon  the 

diffusion  rate  ol  vacancies  T  lie  c  r  i  *  volume,  i  .can 
be  shown  to  depend  on  !  -  in ■<  the  diivmg  lorce 

for  nu.ication  which  is  ,  -  />i.  and  the  site  ol  nu¬ 

cleation  in  the  following  wav 

V  'I 

I  - 1  "  )  1 1 


the  site  of  nucleation.  Note  that  i  is  dimensionless.  It 
has  been  shown  that  the  smallest  values  of  /',  should 
be  achieved  at  junctions  of  particles  and  grain  boun¬ 
daries  [ID],  as  shown  in  Fig.  5  F,  would  depend  on 
the  dihedral  angles  x  a  .  /(  and  /f  all  of  which  must 
satisfy  equilibrium,  while  at  the  same  time,  the  radius 
of  curvature  r  is  uniform  at  all  surfaces  of  the  bound¬ 
ary  The  particle  shape  and  the  surface  energies  will 
influence  F,  which  in  certain  instances  may  even 
become  zero  (lor  example  in  a  matrix  matrix  matrix 
triple  grain  tunction  f,  becomes  zero  if  t  <  r.  bi 
Also.  /,  may  well  vary  from  one  particle  to  another 
giving  nee  to  a  distribution  in  the  values  lor  nu¬ 
cleation  tune  I  xpcriviivc  suggests  that  incubation 
time  would  be  within  the  realm  of  measurement  at 
temperatures  near  o  5  o'  im  when  :  is  ol  the  order  ol 
lie  I o "  atomic  volumes  |  1 1 1  j 

I  he  oxidation  reaction  can  hasten  the  incubation 
time  because  the  cruic.il  nucleus  si/e  will  become 
smwliet  icq.t.iii- •[■  '. <•!  .md  because  the  i.uee  <  O  mol 
codes  m. iv  accclei  ne  cluster  growth  bv  ’.lapping 
vacancies 

4  2  Dtttu  -a  ai  • ,  .j,  r.  ii 

Nucleation  and  growth  ol  yjr.  ;i,o  Kv  involves 

several  sequential  le.tctions  for  ■.•x.imple  as  shown  in 
the  schematic  m  lie  b.  dissociation,  ol  oxygen  to 
atomic  O  at  the  surface  dilhi'ion  o|  (>  alone  "u 
houndarv  to  the  carbide  pre.  .pit  lies  pr  >  idiic  lion  a 
CO  by  oxidation  and  cluster  me  ol  (  t)  wuh  each 
other  and  with  other  vacancies  •  •  t.um  c.ivitie-  and 
diffusion  of  Cl)  back  out  to  the  surtacC  In  principle 
any  one  of  the  above  reactions  can  becon,  late  limit¬ 
ing  Different  scenarios  and  cond, terns  are  surmized 
in  Table  5  The  comments  are  iell  qualitative  became 
at  this  point,  rigorous  analvsis  including  ditlusnm. 
interface  reaction  kinetics  and  cluster  formation 
would  not  lead  to  furthci  insights  Instead  we  wish  to 
emphasize  the  need  tor  more  definitive  and  simple 
experiments 

Some  general  statements  can  be  made  on  the  basis 
of  the  listing  m  I  able  '  High  O  solubility  coupled 
with  low  CO  solubility  in  the  interface  will  promote 
nucleation  and  vice-versa  A  small  carbide  spacing 
itlns  should  be  about  equal  to  the  distance  of  the 


where  is  the  surface  energy.  U„  is  the  weighted  mol¬ 
ecular'*'  volume  and  /  ,  is  a  factor  which  depends  on 

+  Since  the  cacily  mav  contain  vacancies  (  (>  and  GO. 
molecules.  11.  is  the  average  volume  per  molecule  lor  a 
given  composition  of  the  gas' 


fie  '  flic  parameter  / ,  depends  on  the  site  of  nucleation 
ll  is  t lie  smallest  limply  mg  easier  nucleation)  al  the  triple 
junction  between  a  gram  boundary  and  a  second  phase 
particle  ll  depends  on  the  equilibrium  di-hedral  angles  x 
i .  /(  and  />'  when  all  cavity  surfaces  have  uniform  radius  of 
curvature  in  Panicle  shape  and  interface  energies.  Ihcrc- 
lore.  will  influence  F, 
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Fig  6.  Illustration  of  various  kinetic  processes  which  can 
limit  cavity  nucleation. 

carbides  from  the  surface)  will  also  promote  nu¬ 
cleation,  while  presence  of  trace  impurit'es  which  in¬ 
hibit  O  diffusion  will  suppresss  nuclea,ion.  Because 
CO  molecules  are  large  they  would  induce  consider¬ 
able  strain  energy  which  should  relieve  by  having 
vacancies  or  other  CO  molecules  next  to  it.  This 
would  give  rise  to  a  large  binding  energy  and  help 
nucleation.  The  0=0  bond  is  quite  strong,  4.7  eV  and 
there  may  be  instances  where  dissociation  of 
02  — *  20  may  become  important. 

The  solubility  of  O  and  CO  are  likely  to  be 
strongly  influenced  by  the  hydrostatic  component  of 
the  applied  stress  (because  the  partial  molar  volume 
of  solution  may  be  large),  which  means  that  the  stress- 
state  can  affect  environmental  attack.  In  general  a 
tensile  mean  stress,  as  found  at  notch  roots,  will  pro¬ 
mote  nucleation. 

Theoretically,  the  highest  gas  pressures  which  can 
be  achieved  will  be  those  calculated  from  thermodyn¬ 
amics  in  the  previous  section.  In  practice  they  may  be 
lower  because  of  the  kinetic  reasons  described  above. 
Another  limit  on  the  maximum  pressure  is  the  yield 
strength  of  the  material;  if  it  becomes  greater  than  a 
critical  value  then  the  void  will  begin  to  grow  plasti¬ 
cally.  An  approximate  value  for  this  critical  pressure 
can  be  obtained  by  drawing  an  analogy  with  a  thick 
walled  sphere  which  yields  up  to  a  diameter  c  when 
the  cavity  of  diameter,  a,  is  pressurized  to  pressure. 


pmjx.  If  the  material  is  non-hardening  with  a  flow 
stress  V,  then  it  may  be  shown  that  [  2  ^  J 


The  ‘burst’  pressure  can  be  estimated  by  setting  ,  a 
equal  to  (particle  spacing)  (particle  si/ei  Assuming 
the  latter  to  be  about  5.  we  get  that  p,  i  4) 

A  good  discussion  of  the  kinetics  of  void  growth  in 
the  H;  attack  (CH4)  problem  in  steels  can  be  found  in 
Ref.  [24], 


5.  SUMMARY 

Oxidation  of  carbides  (or  carbon)  by  oxygen  trans¬ 
ported  from  the  surface  via  grain  boundaries  can 
produce  CO  at  high  pressures.  This  pressure  adds  to 
tensile  applied  stress  in  promoting  nucleation  and 
growth  of  cavities.  Reasonable  thermodynamic  calcu¬ 
lations  have  been  done,  which  show  that  the  gas 
pressure  can  be  as  large  as  1000  MPa  when  Cr23Cfc 
is  oxidized  at  1000  K  and  pQj  =  1-  At  p0;  =  10“8 
the  pressure  reduces  to  about  100  MPa  and  at 
p0j  =  10" 12  to  about  1  MPa.  From  thermodynamic 
considerations  it  would  appear  that  environmental 
effect  would  begin  at  about  p0j  >  10"  10  bars.  The 
only  data  of  this  type  in  the  literature,  by  Shahinian  et 
al.  [2.  3],  is  in  agreement  with  this.  They  found  that 
low  cycle  fatigue  life  of  Ni  at  823  K  [2]  changes  sig¬ 
moidally  with  increasing  p0j.  The  fatigue  life  first 
begins  to  drop  when  p0,  >  10" 9  and  then  saturates 
when  po,  >  10“ 4  This  would  also  be  expected  from 
the  thermodynamic  calculations:  pco  would  be  much 
smaller  than  the  applied  stress  when  p0,  is  low,  while 
Pco  becomes  large  but  insensitive  to  p0,  when 
Po;  >  10" 4  (see  Fig.  4),  because  fugacity  begins  to  rise 
non-linearly  with  pressure.  The  need  for  further  care¬ 
ful  and  systematic  experimental  work  is  obvious. 

Whereas  the  thermodynamic  criterion  may  provide 
a  qualitative  explanation  for  environmental  attack  in 


Table  5.  Kinetic  limitations  to  cavity  nucleation 


Condition 

Possible  causes 

Slow  0  diffusion 

Trace  impurities  and  traps  in  the  boundary  suppressing  0  diffusion 

y^ 

ii. 

Low  0  solubility 

iii. 

Large  /.  (see  Fig.  6) 

Slow  (interface  limited)  Oxidation  reaction 

Before  nucleation  of  a  cavity 

77 

i. 

Coherent  or  semi-coherent  carbide  matrix  interface 

\rV. 

ii. 

Low  mobility  of  CO  in  the  interface 

After  nucleation  of  a  cavity 

i. 

Structure  of  the  carbide  free  surface 

y 

Sift 

Fa .st  cluster  growth  for  cavity  nucleation 

i. 

High  binding  energy  between  CO  molecules  with  each  other  and 
with  vacancies 

MKN 

ii. 

Low  CO  solubility  in  the  interface 

Slow  CO  diffusion 

i. 

Low  CO  solubility 

ii. 

Trace  impurities  and  traps  in  the  boundary  which  suppress  CO 
diffusion 

iii. 

Large  /  (see  Fig.  6) 
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Ni.  (he  oxygen  attack  in  nickcl-basc  structural  alloys 
appears  to  be  vers  sensitive  to  alloy  composition,  sug¬ 
gesting  the  importance  of  kinetics  At  this  point  it  is 
not  even  possible  to  estimate  whether  nor  not  O  dif¬ 
fusion  is  the  boundarv  is  a  ratc  limiting  step  for  lack 
of  grain  boundarv  diffusion  data  for  O  in  Ni  There  is 
circumstantial  evidence  that  this  may  be  important 
because  addition  of  trace  B  appears  to  hinder  O  dif¬ 
fusion.  thereby  delaying  cavitation  [6] 

7.  ADDITIONAL  REMARKS 

Structural  nickel  base  alloys  contain  a  considerable 
amount  of  C'r  and  Al.  partly  to  enhance  creep  resist¬ 
ance.  but  also  to  provide  a  protective  oxide  on  the 
surface  against  corrosion.  A  tenacious  oxide  layer 
would  reduce  the  oxygen  activity  at  the  oxide-metal 
interface  to  a  near  equilibrium  level,  for  example, 
about  10  30  to  10  21  at  1000  to  I300K  assuming  a 
Cr  Cr2Oj  interface.  Such  low  pa,  would  preclude  the 
cavitation  mechanism  discussed  so  far.  Why  then  are 
structural  alloys  environment  sensitive?  The  question 
can  be  answered  in  several  ways. 

Perhaps  the  most  important  reason  is  that  defor¬ 
mation  can  break  the  surface  layer  permitting  subsur¬ 
face  diffusion  of  Q.  Grain  boundary  sliding  can  pro¬ 
duce  offsets  at  the  surface.  Propagation  of  creep 
cracks,  or  crack-tip  sharpening  in  fatigue  loading,  can 
continually  expose  fresh  metal  surface.  Circumstantial 
support  for  this  idea  comes  from  the  experience  [26] 
that  often  notch  sensitivity  and  environment  sensi¬ 
tivity  of  alloys  goes  hand  in  hand.  Also  alloys  which 
show  stress-rupture  strengthening  in  oxygen  with 
smooth  tensile  specimens,  may  exhibit  oxygen  attack 
if  notched  specimens  are  used.  Presumably  localized 
slip  at  notch  root  breaks  the  oxide  scale.  Experience 
also  suggests  [26]  that  nickel-base  alloys  which  have 
excellent  resistance  to  simple  corrosion,  can  be  quite 
susceptible  to  stress  corrosion  cracking  at  elevated 
temperature. 

In  addition  to  the  internal  gas  production,  internal 
oxidation  of  impurities  or  carbides,  which  produces 
new  oxide  precipitates  or  an  oxide  shell  on  existing 
precipitates  can  also  promote  cavity  nucleation.  The 
reason  is  that,  nearly  always,  oxide  dispersion  alloys 
have  poor  resistance  to  cavitation.  Apparently,  at 
least  in  alloys  of  Cu  and  Ni,  oxide  precipitates  are 
poorly  adhered  to  the  metal.  Internal  oxidation  of 
impurities  in  the  boundaries  could  occur  at  fairly  low 
activities  of  oxygen,  although,  again,  if  a  very  stable 
oxide  such  as  Al203  or  Cr203  forms  on  the  surface 
then  some  mechanism  must  be  available  for  penetrat¬ 
ing  the  oxide  layer. 

There  are  two  other  interesting  observations  in  this 
context.  Oxide  dispersion  strengthened  Ni  alloys 
(Ni-Cr-Th02),  with  equiaxed  grain  structure  are 
known  to  be  much  less  environment  or  notch  sensi¬ 
tive  than  the  precipitation  strengthened  alloys;  pre¬ 
sumably  these  arc  already  prone  to  cavitation  because 
they  contain  oxide  particles  so  that  oxygen  attack  can 


do  no  further  harm  The  second  observation  is  that 
whenever  environmental  attack  is  severe,  the  addition 
of  trace  B  and  / r  has  a  dramatic  beneficial  effect  In 
Ni  C  alloys  this  effect  has  been  correlated  to  the 
obstruction  of  O  diffusion  in  the  presence  of  B;  is  it 
not  possible  that  the  same  mechanism  operates  in  the 
structural  alloys  too'’  I  am  also  led  to  think  that  the 
beneficial  effect  of  B  in  austentic  stainless  steels  which 
cor  tain  some  carbon,  may  also  be  of  a  similar  kind 
.  inally  the  kinetics  of  oxidation  of  metal  surfaces  is 
a  complex  process  even  in  simple  cases  [27.  28],  Our 
current  understanding  of  the  oxidation  of  a  binary 
alloy  such  as  Ni  Cr  is  so  inadequate  that  we  cannot 
even  begin  to  make  judgments  about  the  oxidation 
kinetics  of  muhicomponent  engineering  alloys  [29], 
For  example,  a  Cr,03  layer  on  chromium  may  be 
tenacious  but  in  a  Ni-Cr  alloy  the  oxide  scale  can 
contain  two  or  more  layers  of  different  compositions 
(because  of  kinetic  considerations  [27])  and  therefore, 
it  may  be  more  susceptible  to  cracking  due  to  lattice 
mismatch,  thermal  stresses  and  deformation. 
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